ParaMonte Fortran 2.0.0
Parallel Monte Carlo and Machine Learning Library
See the latest version documentation.
pm_fftnr::getFFTR Interface Reference

Generate and return the Reverse (unnormalized) Fourier Transform of a periodic sequence of type complex or real of arbitrary kind parameter. More...

Detailed Description

Generate and return the Reverse (unnormalized) Fourier Transform of a periodic sequence of type complex or real of arbitrary kind parameter.

See the documentation of pm_fftnr for more details.

Parameters
[in]data: The input contiguous vector of arbitrary size of,
  1. type complex of kind any supported by the processor (e.g., CK, CK32, CK64, or CK128),
  2. type real of kind any supported by the processor (e.g., RK, RK32, RK64, or RK128),
containing the periodic sequence whose FFT is to be computed.
Returns
fft : The output vector of size getExpNext(size(data), 2) of the same type and kind as the input data, containing the FFT result.


Possible calling interfaces

use pm_fftnr, only: getFFTR
type_of(data) :: fft(1 : getExpNext(size(data)))
fft(:) = getFFTR(data(:))
Generate and return the Reverse (unnormalized) Fourier Transform of a periodic sequence of type compl...
Definition: pm_fftnr.F90:611
This module contains procedures and generic interfaces for computing the Discrete Fourier Transform o...
Definition: pm_fftnr.F90:184
Warning
The condition 1 < size(data) must hold for the corresponding input arguments.
This condition is verified only if the library is built with the preprocessor macro CHECK_ENABLED=1.
Remarks
The procedures under discussion are impure.
Note
This functional generic interface is simply a more flexible but slower wrapper around the subroutine generic interface setFFTR.
As such, this functional interface can be significantly slower than the corresponding subroutine interface.
See also
getFFTF
getFFTR
getFFTI
setFFTF
setFFTR
setFFTI


Example usage

1program example
2
3 use pm_kind, only: SK, IK, LK
4 use pm_io, only: display_type
5 use pm_fftnr, only: getFFTF, getFFTR
6 use pm_distUnif, only: getUnifRand
7 use pm_mathCompare, only: isClose
8 use pm_err, only: setAsserted
9
10 implicit none
11
12 integer(IK) :: i
13 integer(IK) :: lenData
14 type(display_type) :: disp
15 disp = display_type(file = "main.out.F90")
16
17 do i = 1, 5
18 block
19 use pm_kind, only: TKC => CKS
20 real(TKC) :: reltol
21 complex(TKC), allocatable :: data(:), atad(:)
22 call disp%show("lenData = getUnifRand(5, 100)")
23 lenData = getUnifRand(5, 100)
24 call disp%show("lenData")
25 call disp%show( lenData )
26 call disp%show("data = 1._TKC + getUnifRand((0._TKC, 0._TKC), (1._TKC, 1._TKC), lenData)")
27 data = 1._TKC + getUnifRand((0._TKC, 0._TKC), (1._TKC, 1._TKC), lenData)
28 call disp%show("data")
29 call disp%show( data )
30 call disp%skip()
31 call disp%show("atad = getFFTR(getFFTF(data))")
32 atad = getFFTR(getFFTF(data))
33 call disp%show("atad = atad / size(atad)")
34 atad = atad / size(atad)
35 call disp%show("atad")
36 call disp%show( atad )
37 call disp%show("reltol = sqrt(epsilon(1._TKC))")
38 reltol = sqrt(epsilon(1._TKC))
39 call disp%show("reltol")
40 call disp%show( reltol )
41 call disp%show("isClose(data, atad(1:size(data)), reltol = reltol)")
42 call disp%show( isClose(data, atad(1:size(data)), reltol = reltol) )
43 call disp%show("call setAsserted(all(isClose(data, atad(1:size(data)), reltol = reltol)))")
44 call setAsserted(all(isClose(data, atad(1:size(data)), reltol = reltol)))
45 call disp%skip()
46 end block
47
48 block
49 use pm_kind, only: TKC => CKD
50 real(TKC) :: reltol
51 complex(TKC), allocatable :: data(:), atad(:)
52 call disp%show("lenData = getUnifRand(5, 100)")
53 lenData = getUnifRand(5, 100)
54 call disp%show("lenData")
55 call disp%show( lenData )
56 call disp%show("data = 1._TKC + getUnifRand((0._TKC, 0._TKC), (1._TKC, 1._TKC), lenData)")
57 data = 1._TKC + getUnifRand((0._TKC, 0._TKC), (1._TKC, 1._TKC), lenData)
58 call disp%show("data")
59 call disp%show( data )
60 call disp%skip()
61 call disp%show("atad = getFFTR(getFFTF(data))")
62 atad = getFFTR(getFFTF(data))
63 call disp%show("atad = atad / size(atad)")
64 atad = atad / size(atad)
65 call disp%show("atad")
66 call disp%show( atad )
67 call disp%show("reltol = sqrt(epsilon(1._TKC))")
68 reltol = sqrt(epsilon(1._TKC))
69 call disp%show("reltol")
70 call disp%show( reltol )
71 call disp%show("isClose(data, atad(1:size(data)), reltol = reltol)")
72 call disp%show( isClose(data, atad(1:size(data)), reltol = reltol) )
73 call disp%show("call setAsserted(all(isClose(data, atad(1:size(data)), reltol = reltol)))")
74 call setAsserted(all(isClose(data, atad(1:size(data)), reltol = reltol)))
75 call disp%skip()
76 end block
77
78 block
79 use pm_kind, only: TKC => CKH
80 real(TKC) :: reltol
81 complex(TKC), allocatable :: data(:), atad(:)
82 call disp%show("lenData = getUnifRand(5, 100)")
83 lenData = getUnifRand(5, 100)
84 call disp%show("lenData")
85 call disp%show( lenData )
86 call disp%show("data = 1._TKC + getUnifRand((0._TKC, 0._TKC), (1._TKC, 1._TKC), lenData)")
87 data = 1._TKC + getUnifRand((0._TKC, 0._TKC), (1._TKC, 1._TKC), lenData)
88 call disp%show("data")
89 call disp%show( data )
90 call disp%skip()
91 call disp%show("atad = getFFTR(getFFTF(data))")
92 atad = getFFTR(getFFTF(data))
93 call disp%show("atad = atad / size(atad)")
94 atad = atad / size(atad)
95 call disp%show("atad")
96 call disp%show( atad )
97 call disp%show("reltol = sqrt(epsilon(1._TKC))")
98 reltol = sqrt(epsilon(1._TKC))
99 call disp%show("reltol")
100 call disp%show( reltol )
101 call disp%show("isClose(data, atad(1:size(data)), reltol = reltol)")
102 call disp%show( isClose(data, atad(1:size(data)), reltol = reltol) )
103 call disp%show("call setAsserted(all(isClose(data, atad(1:size(data)), reltol = reltol)))")
104 call setAsserted(all(isClose(data, atad(1:size(data)), reltol = reltol)))
105 call disp%skip()
106 end block
107
108 block
109 use pm_kind, only: TKC => RKS
110 real(TKC) :: reltol
111 real(TKC), allocatable :: data(:), atad(:)
112 call disp%show("lenData = getUnifRand(5, 100)")
113 lenData = getUnifRand(5, 100)
114 call disp%show("lenData")
115 call disp%show( lenData )
116 call disp%show("data = 1._TKC + getUnifRand(0._TKC, 1._TKC, lenData)")
117 data = 1._TKC + getUnifRand(0._TKC, 1._TKC, lenData)
118 call disp%show("data")
119 call disp%show( data )
120 call disp%skip()
121 call disp%show("atad = getFFTR(getFFTF(data))")
122 atad = getFFTR(getFFTF(data))
123 call disp%show("atad = atad * 2 / size(atad)")
124 atad = atad * 2 / size(atad)
125 call disp%show("atad")
126 call disp%show( atad )
127 call disp%show("reltol = sqrt(epsilon(1._TKC))")
128 reltol = sqrt(epsilon(1._TKC))
129 call disp%show("reltol")
130 call disp%show( reltol )
131 call disp%show("isClose(data, atad(1:size(data)), reltol = reltol)")
132 call disp%show( isClose(data, atad(1:size(data)), reltol = reltol) )
133 call disp%show("call setAsserted(all(isClose(data, atad(1:size(data)), reltol = reltol)))")
134 call setAsserted(all(isClose(data, atad(1:size(data)), reltol = reltol)))
135 call disp%skip()
136 end block
137
138 block
139 use pm_kind, only: TKC => RKD
140 real(TKC) :: reltol
141 real(TKC), allocatable :: data(:), atad(:)
142 call disp%show("lenData = getUnifRand(5, 100)")
143 lenData = getUnifRand(5, 100)
144 call disp%show("lenData")
145 call disp%show( lenData )
146 call disp%show("data = 1._TKC + getUnifRand(0._TKC, 1._TKC, lenData)")
147 data = 1._TKC + getUnifRand(0._TKC, 1._TKC, lenData)
148 call disp%show("data")
149 call disp%show( data )
150 call disp%skip()
151 call disp%show("atad = getFFTR(getFFTF(data))")
152 atad = getFFTR(getFFTF(data))
153 call disp%show("atad = atad * 2 / size(atad)")
154 atad = atad * 2 / size(atad)
155 call disp%show("atad")
156 call disp%show( atad )
157 call disp%show("reltol = sqrt(epsilon(1._TKC))")
158 reltol = sqrt(epsilon(1._TKC))
159 call disp%show("reltol")
160 call disp%show( reltol )
161 call disp%show("isClose(data, atad(1:size(data)), reltol = reltol)")
162 call disp%show( isClose(data, atad(1:size(data)), reltol = reltol) )
163 call disp%show("call setAsserted(all(isClose(data, atad(1:size(data)), reltol = reltol)))")
164 call setAsserted(all(isClose(data, atad(1:size(data)), reltol = reltol)))
165 call disp%skip()
166 end block
167
168 block
169 use pm_kind, only: TKC => RKH
170 real(TKC) :: reltol
171 real(TKC), allocatable :: data(:), atad(:)
172 call disp%show("lenData = getUnifRand(5, 100)")
173 lenData = getUnifRand(5, 100)
174 call disp%show("lenData")
175 call disp%show( lenData )
176 call disp%show("data = 1._TKC + getUnifRand(0._TKC, 1._TKC, lenData)")
177 data = 1._TKC + getUnifRand(0._TKC, 1._TKC, lenData)
178 call disp%show("data")
179 call disp%show( data )
180 call disp%skip()
181 call disp%show("atad = getFFTR(getFFTF(data))")
182 atad = getFFTR(getFFTF(data))
183 call disp%show("atad = atad * 2 / size(atad)")
184 atad = atad * 2 / size(atad)
185 call disp%show("atad")
186 call disp%show( atad )
187 call disp%show("reltol = sqrt(epsilon(1._TKC))")
188 reltol = sqrt(epsilon(1._TKC))
189 call disp%show("reltol")
190 call disp%show( reltol )
191 call disp%show("isClose(data, atad(1:size(data)), reltol = reltol)")
192 call disp%show( isClose(data, atad(1:size(data)), reltol = reltol) )
193 call disp%show("call setAsserted(all(isClose(data, atad(1:size(data)), reltol = reltol)))")
194 call setAsserted(all(isClose(data, atad(1:size(data)), reltol = reltol)))
195 call disp%skip()
196 end block
197 end do
198
199end program example
Generate and return a scalar or a contiguous array of rank 1 of length s1 of randomly uniformly distr...
Verify the input assertion holds and if it does not, print the (optional) input message on stdout and...
Definition: pm_err.F90:735
Generate and return an object of type stop_type with the user-specified input attributes.
Definition: pm_err.F90:1618
Generate and return the Forward Fourier Transform (a.k.a. Fourier Analysis) of a periodic sequence of...
Definition: pm_fftnr.F90:250
This is a generic method of the derived type display_type with pass attribute.
Definition: pm_io.F90:11131
Generate and return .true. if the two input values are sufficiently close to each other within the sp...
This module contains classes and procedures for computing various statistical quantities related to t...
This module contains classes and procedures for reporting and handling errors.
Definition: pm_err.F90:52
This module contains classes and procedures for input/output (IO) or generic display operations on st...
Definition: pm_io.F90:252
type(display_type) disp
This is a scalar module variable an object of type display_type for general display.
Definition: pm_io.F90:11016
This module defines the relevant Fortran kind type-parameters frequently used in the ParaMonte librar...
Definition: pm_kind.F90:268
integer, parameter LK
The default logical kind in the ParaMonte library: kind(.true.) in Fortran, kind(....
Definition: pm_kind.F90:541
integer, parameter CKH
The scalar integer constant of intrinsic default kind, representing the highest-precision complex kin...
Definition: pm_kind.F90:843
integer, parameter CKS
The single-precision complex kind in Fortran mode. On most platforms, this is a 32-bit real kind.
Definition: pm_kind.F90:570
integer, parameter IK
The default integer kind in the ParaMonte library: int32 in Fortran, c_int32_t in C-Fortran Interoper...
Definition: pm_kind.F90:540
integer, parameter CKD
The double precision complex kind in Fortran mode. On most platforms, this is a 64-bit real kind.
Definition: pm_kind.F90:571
integer, parameter RKD
The double precision real kind in Fortran mode. On most platforms, this is an 64-bit real kind.
Definition: pm_kind.F90:568
integer, parameter SK
The default character kind in the ParaMonte library: kind("a") in Fortran, c_char in C-Fortran Intero...
Definition: pm_kind.F90:539
integer, parameter RKH
The scalar integer constant of intrinsic default kind, representing the highest-precision real kind t...
Definition: pm_kind.F90:858
integer, parameter RKS
The single-precision real kind in Fortran mode. On most platforms, this is an 32-bit real kind.
Definition: pm_kind.F90:567
This module contains the procedures and interfaces for evaluating the relative or absolute proximity ...
Generate and return an object of type display_type.
Definition: pm_io.F90:10282

Example Unix compile command via Intel ifort compiler
1#!/usr/bin/env sh
2rm main.exe
3ifort -fpp -standard-semantics -O3 -Wl,-rpath,../../../lib -I../../../inc main.F90 ../../../lib/libparamonte_fortran_*_intel* -o main.exe
4./main.exe

Example Windows Batch compile command via Intel ifort compiler
1del main.exe
2set PATH=..\..\..\lib;%PATH%
3ifort /fpp /standard-semantics /O3 /I:..\..\..\include main.F90 ..\..\..\lib\libparamonte_fortran_*_intel*.lib /exe:main.exe
4main.exe

Example Unix / MinGW compile command via GNU gfortran compiler
1#!/usr/bin/env sh
2rm main.exe
3gfortran -cpp -ffree-line-length-none -O3 -Wl,-rpath,../../../lib -I../../../inc main.F90 ../../../lib/libparamonte_fortran_*_gnu* -o main.exe
4./main.exe

Example output
1lenData = getUnifRand(5, 100)
2lenData
3+39
4data = 1._TKC + getUnifRand((0._TKC, 0._TKC), (1._TKC, 1._TKC), lenData)
5data
6(+1.42792273, +0.181441665), (+1.28003955, +0.739739835), (+1.48283458, +0.395718217), (+1.58713222, +0.433815181), (+1.81654048, +0.735282660), (+1.71164989, +0.897854507), (+1.97664309, +0.419615090), (+1.08178687, +0.177605212), (+1.92551184, +0.394518256), (+1.39148295, +0.942302942), (+1.37600565, +0.196714938), (+1.15786839, +0.971206069), (+1.97886467, +0.546871722), (+1.75491512, +0.605235994), (+1.12047029, +0.873045921), (+1.16300225, +0.982885420), (+1.84586525, +0.121501505), (+1.68181467, +0.559351683), (+1.75846076, +0.775392473), (+1.48434091, +0.101965547), (+1.78746867, +0.544283628), (+1.98398733, +0.620429218), (+1.50212812, +0.325695634), (+1.02254367, +0.863959908), (+1.20206714, +0.210701644), (+1.47501862, +0.680518627), (+1.29716253, +0.873645842), (+1.84417009, +0.972444177), (+1.18429375, +0.119709611), (+1.68508673, +0.287689567), (+1.10278440, +0.322481632), (+1.62978697, +0.832166910), (+1.23114920, +0.214612424), (+1.22124910, +0.251632333), (+1.19969606, +0.314506292E-1), (+1.02864099, +0.108356595), (+1.99477136, +0.161579669), (+1.59786725, +0.729205906), (+1.97374845, +0.736678958)
7
8atad = getFFTR(getFFTF(data))
9atad = atad / size(atad)
10atad
11(+1.42792273, +0.181441709), (+1.28003967, +0.739739895), (+1.48283434, +0.395718217), (+1.58713186, +0.433814943), (+1.81654048, +0.735282779), (+1.71164978, +0.897854567), (+1.97664320, +0.419615120), (+1.08178675, +0.177605078), (+1.92551160, +0.394518435), (+1.39148259, +0.942303061), (+1.37600565, +0.196714953), (+1.15786827, +0.971205950), (+1.97886443, +0.546872020), (+1.75491500, +0.605236173), (+1.12047017, +0.873046100), (+1.16300225, +0.982885361), (+1.84586513, +0.121501431), (+1.68181455, +0.559351563), (+1.75846064, +0.775392234), (+1.48434043, +0.101965442), (+1.78746867, +0.544283569), (+1.98398721, +0.620429218), (+1.50212812, +0.325695574), (+1.02254343, +0.863959849), (+1.20206690, +0.210701555), (+1.47501838, +0.680518508), (+1.29716253, +0.873645663), (+1.84416986, +0.972443938), (+1.18429363, +0.119709559), (+1.68508649, +0.287689537), (+1.10278440, +0.322481573), (+1.62978697, +0.832166910), (+1.23114920, +0.214612469), (+1.22124922, +0.251632214), (+1.19969606, +0.314505994E-1), (+1.02864110, +0.108356670), (+1.99477124, +0.161579728), (+1.59786713, +0.729205847), (+1.97374856, +0.736678839), (+0.119209290E-6, +0.223517418E-6), (+0.119209290E-6, -0.745058060E-7), (+0.298023224E-6, -0.894069672E-7), (+0.417232513E-6, -0.149011612E-7), (+0.119209290E-6, +0.298023224E-7), (+0.596046448E-7, -0.178813934E-6), (+0.00000000, -0.178813934E-6), (+0.119209290E-6, +0.596046448E-7), (+0.596046448E-7, -0.298023224E-7), (+0.119209290E-6, +0.447034836E-7), (+0.238418579E-6, +0.894069672E-7), (+0.238418579E-6, +0.178813934E-6), (-0.596046448E-7, +0.193715096E-6), (+0.119209290E-6, +0.596046448E-7), (+0.119209290E-6, +0.596046448E-7), (+0.119209290E-6, +0.447034836E-7), (+0.00000000, +0.178813934E-6), (+0.178813934E-6, +0.119209290E-6), (+0.238418579E-6, +0.894069672E-7), (+0.357627869E-6, +0.149011612E-6), (+0.238418579E-6, +0.208616257E-6), (+0.238418579E-6, +0.372529030E-7), (+0.298023224E-6, +0.298023224E-7), (+0.596046448E-7, +0.745058060E-7), (+0.00000000, +0.149011612E-6)
12reltol = sqrt(epsilon(1._TKC))
13reltol
14+0.345266977E-3
15isClose(data, atad(1:size(data)), reltol = reltol)
16T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T
17call setAsserted(all(isClose(data, atad(1:size(data)), reltol = reltol)))
18
19lenData = getUnifRand(5, 100)
20lenData
21+38
22data = 1._TKC + getUnifRand((0._TKC, 0._TKC), (1._TKC, 1._TKC), lenData)
23data
24(+1.1937012649361112, +0.84427392435013393), (+1.1953154210883299, +0.76399798932981455), (+1.4781020575041377, +0.50347765606598993), (+1.0677930629617935, +0.62296375636851575), (+1.2210525337665312, +0.41049907221561044), (+1.1207037468822814, +0.13732587815357133), (+1.1001783349986485, +0.99930371283357755), (+1.8442553612035173, +0.99230586044953129), (+1.7004451212962461, +0.45962005578628085), (+1.5320954448580335, +0.31980527085163080E-1), (+1.6090440514861690, +0.79093128586185069), (+1.0742767429820348, +0.71050999235403889), (+1.6503180980175778, +0.49312299142385629), (+1.7786641045966189, +0.62218001553465263), (+1.8993772476683910, +0.84167476703835375), (+1.5404822198849830, +0.12793493979111303), (+1.5032681613897843, +0.66756464940626004), (+1.7113958656642438, +0.75268880692084006), (+1.9494086257104688, +0.32356565689564654), (+1.9186817246872367, +0.53575284702563009E-1), (+1.6898670190862859, +0.83894965448797032), (+1.5866897610842154, +0.44715265569456408), (+1.4992226558432082, +0.24304598562694357), (+1.8025908775314208, +0.49886374654192844), (+1.3361560698894765, +0.89125762607055214), (+1.0051069298470514, +0.60874704540627478), (+1.4551851252705474, +0.73544623672448717), (+1.0254675599234506, +0.84263150714910162), (+1.3935307634780951, +0.82424983406904295E-1), (+1.8762961774542726, +0.65366351258824540), (+1.3465804901516891, +0.42463493595189838), (+1.1416094092706657, +0.45859094966813652), (+1.6348319291216613, +0.71226990878158647E-1), (+1.1986367163334211, +0.73855240655769316), (+1.4157286404415732, +0.87380544524790504), (+1.0495194701515300, +0.48528619103817516), (+1.3506870715328794, +0.59251991521067560), (+1.2086518974047544, +0.73704525070437887)
25
26atad = getFFTR(getFFTF(data))
27atad = atad / size(atad)
28atad
29(+1.1937012649361112, +0.84427392435013382), (+1.1953154210883299, +0.76399798932981466), (+1.4781020575041375, +0.50347765606599015), (+1.0677930629617931, +0.62296375636851586), (+1.2210525337665312, +0.41049907221561038), (+1.1207037468822816, +0.13732587815357139), (+1.1001783349986487, +0.99930371283357755), (+1.8442553612035173, +0.99230586044953140), (+1.7004451212962461, +0.45962005578628096), (+1.5320954448580333, +0.31980527085163440E-1), (+1.6090440514861686, +0.79093128586185069), (+1.0742767429820350, +0.71050999235403900), (+1.6503180980175776, +0.49312299142385635), (+1.7786641045966189, +0.62218001553465285), (+1.8993772476683910, +0.84167476703835375), (+1.5404822198849828, +0.12793493979111362), (+1.5032681613897843, +0.66756464940625992), (+1.7113958656642438, +0.75268880692083995), (+1.9494086257104688, +0.32356565689564609), (+1.9186817246872367, +0.53575284702562842E-1), (+1.6898670190862861, +0.83894965448797021), (+1.5866897610842159, +0.44715265569456419), (+1.4992226558432085, +0.24304598562694310), (+1.8025908775314208, +0.49886374654192811), (+1.3361560698894763, +0.89125762607055170), (+1.0051069298470514, +0.60874704540627467), (+1.4551851252705474, +0.73544623672448628), (+1.0254675599234511, +0.84263150714910107), (+1.3935307634780951, +0.82424983406904087E-1), (+1.8762961774542726, +0.65366351258824529), (+1.3465804901516891, +0.42463493595189772), (+1.1416094092706657, +0.45859094966813585), (+1.6348319291216611, +0.71226990878158758E-1), (+1.1986367163334211, +0.73855240655769305), (+1.4157286404415734, +0.87380544524790404), (+1.0495194701515302, +0.48528619103817483), (+1.3506870715328794, +0.59251991521067560), (+1.2086518974047542, +0.73704525070437876), (+0.44408920985006262E-15, -0.27755575615628914E-15), (+0.22204460492503131E-15, -0.38857805861880479E-15), (-0.44408920985006262E-15, -0.19428902930940239E-15), (+0.22204460492503131E-15, -0.30531133177191805E-15), (+0.55511151231257827E-15, -0.55511151231257827E-16), (+0.22204460492503131E-15, -0.11102230246251565E-15), (+0.11102230246251565E-15, -0.27755575615628914E-15), (+0.33306690738754696E-15, -0.38857805861880479E-15), (+0.55511151231257827E-15, -0.16653345369377348E-15), (+0.22204460492503131E-15, -0.47184478546569153E-15), (-0.22204460492503131E-15, +0.16653345369377348E-15), (-0.11102230246251565E-15, +0.16653345369377348E-15), (+0.22204460492503131E-15, +0.58286708792820718E-15), (+0.11102230246251565E-15, +0.55511151231257827E-16), (-0.22204460492503131E-15, +0.55511151231257827E-16), (-0.11102230246251565E-15, -0.11102230246251565E-15), (+0.0000000000000000, +0.52735593669694936E-15), (+0.22204460492503131E-15, +0.22204460492503131E-15), (-0.11102230246251565E-15, +0.11102230246251565E-15), (+0.11102230246251565E-15, +0.0000000000000000), (+0.11102230246251565E-15, +0.83266726846886741E-15), (+0.11102230246251565E-15, +0.38857805861880479E-15), (+0.0000000000000000, +0.13877787807814457E-16), (+0.22204460492503131E-15, +0.55511151231257827E-16), (+0.66613381477509392E-15, +0.99920072216264089E-15), (+0.33306690738754696E-15, +0.55511151231257827E-15)
30reltol = sqrt(epsilon(1._TKC))
31reltol
32+0.14901161193847656E-7
33isClose(data, atad(1:size(data)), reltol = reltol)
34T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T
35call setAsserted(all(isClose(data, atad(1:size(data)), reltol = reltol)))
36
37lenData = getUnifRand(5, 100)
38lenData
39+38
40data = 1._TKC + getUnifRand((0._TKC, 0._TKC), (1._TKC, 1._TKC), lenData)
41data
42(+1.41220134025324635844386276159111801, +0.548322190832395844851194478685682877), (+1.50794109305555591980276828446173283, +0.515705460382486012843437261589866987), (+1.40986980768943715157328374715022141, +0.128051299569145666997101067633448842), (+1.26542549425400164599721673445397319, +0.697446042937342457092477368936398343), (+1.47182430261348232532083550539919805, +0.999437203572338303622136490707201134), (+1.83897147232745391209493204194630956, +0.318713276662000535617713289652920344), (+1.64692726079628702482950655761038190, +0.126971105676949314131754194674744755), (+1.26019721250719591462007826370901776, +0.705486776828076622006212099322418258), (+1.92797430928329383446312238695364366, +0.157342559191610247686434035369926568), (+1.02154058752549187883899604780025786, +0.413688094961745262947574867017725298), (+1.66834898852583973201478841494207976, +0.492841162550939994103755373953696084E-1), (+1.26984857248265775118134621738085328, +0.861569294123515426496483240700626627), (+1.67518681062239268232653000195969670, +0.703514483576264174089176226523405022E-1), (+1.12858444481571292366835353735552589, +0.911547691252690101606472811204564345), (+1.74126198541552277843427213689270066, +0.312042843575847471917055075985806169), (+1.14160014644244172757072165076983265, +0.119883414746902710359271744898532428E-1), (+1.53649584316602925998271035975892402, +0.950591478560959946898079574718849175), (+1.82109033622705772584175918626879078, +0.697328543829916435446722255639953613), (+1.04366181768963701298540426008620741, +0.765927542033013462627522454880685693), (+1.32289721590781302410177221925148077, +0.573113446073010923060297285307612792), (+1.44598445545785295143379956486921021, +0.451390424640126981945888450055901526), (+1.14529899007379108115157344057937685, +0.226151601571136174249685882964649126), (+1.31895499933540964031541596219388963, +0.139976831455383365783668714457847194), (+1.09807295510050796529756201978916781, +0.989488382505176125613450219517983727E-1), (+1.23102067001048198725637241423528414, +0.844363496907115040306763739698050744), (+1.02942619070398029133860540606685333, +0.716743220880416637807721670043283975), (+1.41755275060014360453080117573292708, +0.566846702001763961129338078031366478), (+1.92005596639579703536655856051464133, +0.867765664778959766110789117599135713), (+1.76285613905381861222653116311365506, +0.663984498888541357492765118312060462), (+1.58020293486287711377540248066352866, +0.320698151145023918200650037899173329), (+1.01172400411417005316729129332185964, +0.674658992387062601138748650908601046), (+1.03212486940806763290536279293883070, +0.533445912366626459799927646724086821), (+1.69691817118665690992211559375194824, +0.490599418838769862905634431248622481), (+1.73296086452403526703797622979282048, +0.150123573837736280834063758116625616), (+1.64133553060108381807742152171065401, +0.819859184009136321743377446741342416), (+1.55706354265378620805291684490730865, +0.272462605056405761768117139072994197E-2), (+1.30929650149634923126287939013707574, +0.196635336473957644880202260604521101E-1), (+1.16920054513906050423583528198122210, +0.700126872762604166892415873062181023)
43
44atad = getFFTR(getFFTF(data))
45atad = atad / size(atad)
46atad
47(+1.41220134025324635844386276159111782, +0.548322190832395844851194478685682974), (+1.50794109305555591980276828446173302, +0.515705460382486012843437261589867276), (+1.40986980768943715157328374715022141, +0.128051299569145666997101067633449082), (+1.26542549425400164599721673445397338, +0.697446042937342457092477368936398535), (+1.47182430261348232532083550539919805, +0.999437203572338303622136490707201230), (+1.83897147232745391209493204194631033, +0.318713276662000535617713289652920247), (+1.64692726079628702482950655761038190, +0.126971105676949314131754194674744562), (+1.26019721250719591462007826370901814, +0.705486776828076622006212099322418355), (+1.92797430928329383446312238695364366, +0.157342559191610247686434035369926664), (+1.02154058752549187883899604780025786, +0.413688094961745262947574867017725490), (+1.66834898852583973201478841494207938, +0.492841162550939994103755373953696325E-1), (+1.26984857248265775118134621738085308, +0.861569294123515426496483240700626820), (+1.67518681062239268232653000195969670, +0.703514483576264174089176226523407429E-1), (+1.12858444481571292366835353735552608, +0.911547691252690101606472811204564923), (+1.74126198541552277843427213689270104, +0.312042843575847471917055075985806361), (+1.14160014644244172757072165076983265, +0.119883414746902710359271744898533451E-1), (+1.53649584316602925998271035975892422, +0.950591478560959946898079574718849175), (+1.82109033622705772584175918626879078, +0.697328543829916435446722255639953709), (+1.04366181768963701298540426008620741, +0.765927542033013462627522454880685597), (+1.32289721590781302410177221925148077, +0.573113446073010923060297285307612600), (+1.44598445545785295143379956486921060, +0.451390424640126981945888450055901526), (+1.14529899007379108115157344057937685, +0.226151601571136174249685882964649006), (+1.31895499933540964031541596219388963, +0.139976831455383365783668714457847145), (+1.09807295510050796529756201978916800, +0.989488382505176125613450219517982283E-1), (+1.23102067001048198725637241423528414, +0.844363496907115040306763739698050744), (+1.02942619070398029133860540606685313, +0.716743220880416637807721670043283975), (+1.41755275060014360453080117573292689, +0.566846702001763961129338078031365997), (+1.92005596639579703536655856051464095, +0.867765664778959766110789117599135232), (+1.76285613905381861222653116311365467, +0.663984498888541357492765118312060462), (+1.58020293486287711377540248066352885, +0.320698151145023918200650037899173233), (+1.01172400411417005316729129332185983, +0.674658992387062601138748650908601143), (+1.03212486940806763290536279293883109, +0.533445912366626459799927646724086628), (+1.69691817118665690992211559375194824, +0.490599418838769862905634431248622529), (+1.73296086452403526703797622979282029, +0.150123573837736280834063758116625447), (+1.64133553060108381807742152171065440, +0.819859184009136321743377446741342224), (+1.55706354265378620805291684490730885, +0.272462605056405761768117139072979752E-2), (+1.30929650149634923126287939013707574, +0.196635336473957644880202260604522064E-1), (+1.16920054513906050423583528198122210, +0.700126872762604166892415873062181216), (+0.00000000000000000000000000000000000, +0.00000000000000000000000000000000000), (+0.962964972193617926527988971292463659E-34, -0.288889491658085377958396691387739098E-33), (+0.962964972193617926527988971292463659E-34, +0.00000000000000000000000000000000000), (+0.962964972193617926527988971292463659E-34, -0.216667118743564033468797518540804323E-33), (+0.192592994438723585305597794258492732E-33, -0.168518870133883137142398069976181140E-33), (+0.192592994438723585305597794258492732E-33, -0.144444745829042688979198345693869549E-33), (-0.962964972193617926527988971292463659E-34, -0.240741243048404481631997242823115915E-33), (-0.192592994438723585305597794258492732E-33, -0.385185988877447170611195588516985464E-33), (+0.962964972193617926527988971292463659E-34, -0.240741243048404481631997242823115915E-34), (+0.962964972193617926527988971292463659E-34, -0.150463276905252801019998276764447447E-33), (+0.192592994438723585305597794258492732E-33, +0.00000000000000000000000000000000000), (+0.962964972193617926527988971292463659E-34, -0.962964972193617926527988971292463659E-34), (+0.192592994438723585305597794258492732E-33, +0.962964972193617926527988971292463659E-34), (+0.288889491658085377958396691387739098E-33, +0.962964972193617926527988971292463659E-34), (-0.288889491658085377958396691387739098E-33, -0.722223729145213444895991728469347744E-34), (+0.00000000000000000000000000000000000, -0.722223729145213444895991728469347744E-34), (-0.962964972193617926527988971292463659E-34, +0.144444745829042688979198345693869549E-33), (+0.192592994438723585305597794258492732E-33, +0.00000000000000000000000000000000000), (+0.192592994438723585305597794258492732E-33, +0.00000000000000000000000000000000000), (+0.192592994438723585305597794258492732E-33, +0.962964972193617926527988971292463659E-34), (+0.00000000000000000000000000000000000, +0.385185988877447170611195588516985464E-33), (+0.00000000000000000000000000000000000, +0.240741243048404481631997242823115915E-33), (+0.962964972193617926527988971292463659E-34, -0.962964972193617926527988971292463659E-34), (+0.962964972193617926527988971292463659E-34, -0.288889491658085377958396691387739098E-33), (+0.00000000000000000000000000000000000, +0.00000000000000000000000000000000000), (+0.288889491658085377958396691387739098E-33, +0.00000000000000000000000000000000000)
48reltol = sqrt(epsilon(1._TKC))
49reltol
50+0.138777878078144567552953958511352539E-16
51isClose(data, atad(1:size(data)), reltol = reltol)
52T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T
53call setAsserted(all(isClose(data, atad(1:size(data)), reltol = reltol)))
54
55lenData = getUnifRand(5, 100)
56lenData
57+6
58data = 1._TKC + getUnifRand(0._TKC, 1._TKC, lenData)
59data
60+1.24285269, +1.04180515, +1.16328835, +1.94218338, +1.99697053, +1.92897630
61
62atad = getFFTR(getFFTF(data))
63atad = atad * 2 / size(atad)
64atad
65+1.24285269, +1.04180515, +1.16328824, +1.94218326, +1.99697065, +1.92897618, +0.00000000, +0.178813934E-6
66reltol = sqrt(epsilon(1._TKC))
67reltol
68+0.345266977E-3
69isClose(data, atad(1:size(data)), reltol = reltol)
70T, T, T, T, T, T
71call setAsserted(all(isClose(data, atad(1:size(data)), reltol = reltol)))
72
73lenData = getUnifRand(5, 100)
74lenData
75+46
76data = 1._TKC + getUnifRand(0._TKC, 1._TKC, lenData)
77data
78+1.8403126972916042, +1.6262086221388476, +1.7545219984537272, +1.9615826523911757, +1.5492810800075216, +1.0316834947438527, +1.1888889323015803, +1.7158951436480820, +1.2083960870078692, +1.3102120495802487, +1.4448444269942802, +1.3170423959539308, +1.9970818616088928, +1.8899449898147667, +1.8380226319053390, +1.3868220094295955, +1.2982288477899568, +1.1538849751134508, +1.7449855718227276, +1.9667468555698859, +1.2347801983981086, +1.6337576773060163, +1.5058123929073144, +1.5863231919046594, +1.2731162816641775, +1.1538316568381677, +1.2676645271914313, +1.1389761543392256, +1.1581031182506276, +1.8169295730915389, +1.6895101419265643, +1.4441664151119298, +1.8678555982239886, +1.2263543162184862, +1.5307520213427961, +1.0724334101141544, +1.1568046255217816, +1.5322028390940221, +1.2002480359441443, +1.1102048372257918, +1.9994050256989193, +1.3820783781964645, +1.7593451434829157, +1.9684822739487136, +1.4750138105508137, +1.3301957275432601
79
80atad = getFFTR(getFFTF(data))
81atad = atad * 2 / size(atad)
82atad
83+1.8403126972916040, +1.6262086221388479, +1.7545219984537268, +1.9615826523911757, +1.5492810800075214, +1.0316834947438527, +1.1888889323015799, +1.7158951436480820, +1.2083960870078689, +1.3102120495802487, +1.4448444269942800, +1.3170423959539306, +1.9970818616088923, +1.8899449898147664, +1.8380226319053381, +1.3868220094295953, +1.2982288477899566, +1.1538849751134506, +1.7449855718227274, +1.9667468555698855, +1.2347801983981086, +1.6337576773060163, +1.5058123929073144, +1.5863231919046590, +1.2731162816641772, +1.1538316568381677, +1.2676645271914313, +1.1389761543392249, +1.1581031182506276, +1.8169295730915382, +1.6895101419265641, +1.4441664151119291, +1.8678555982239886, +1.2263543162184860, +1.5307520213427961, +1.0724334101141544, +1.1568046255217814, +1.5322028390940226, +1.2002480359441443, +1.1102048372257918, +1.9994050256989189, +1.3820783781964647, +1.7593451434829157, +1.9684822739487131, +1.4750138105508142, +1.3301957275432601, +0.11102230246251565E-14, -0.11102230246251565E-15, +0.22204460492503131E-15, +0.11102230246251565E-15, +0.44408920985006262E-15, +0.22204460492503131E-15, +0.22204460492503131E-15, +0.11102230246251565E-15, +0.55511151231257827E-15, +0.55511151231257827E-15, +0.11102230246251565E-15, +0.11102230246251565E-15, +0.55511151231257827E-15, +0.77715611723760958E-15, +0.44408920985006262E-15, +0.11102230246251565E-15, +0.33306690738754696E-15, +0.88817841970012523E-15
84reltol = sqrt(epsilon(1._TKC))
85reltol
86+0.14901161193847656E-7
87isClose(data, atad(1:size(data)), reltol = reltol)
88T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T
89call setAsserted(all(isClose(data, atad(1:size(data)), reltol = reltol)))
90
91lenData = getUnifRand(5, 100)
92lenData
93+8
94data = 1._TKC + getUnifRand(0._TKC, 1._TKC, lenData)
95data
96+1.04871433464436812316952739328355120, +1.94471551166988120603155587794268733, +1.43162914763762258904747324322788045, +1.74600098035463228302072437767348585, +1.29619953489561099207123680760600153, +1.37636200874256764945640196114032300, +1.93353367852645537478682051658537684, +1.23992748549319634559587057127249249
97
98atad = getFFTR(getFFTF(data))
99atad = atad * 2 / size(atad)
100atad
101+1.04871433464436812316952739328355100, +1.94471551166988120603155587794268733, +1.43162914763762258904747324322788045, +1.74600098035463228302072437767348585, +1.29619953489561099207123680760600134, +1.37636200874256764945640196114032300, +1.93353367852645537478682051658537665, +1.23992748549319634559587057127249249
102reltol = sqrt(epsilon(1._TKC))
103reltol
104+0.138777878078144567552953958511352539E-16
105isClose(data, atad(1:size(data)), reltol = reltol)
106T, T, T, T, T, T, T, T
107call setAsserted(all(isClose(data, atad(1:size(data)), reltol = reltol)))
108
109lenData = getUnifRand(5, 100)
110lenData
111+93
112data = 1._TKC + getUnifRand((0._TKC, 0._TKC), (1._TKC, 1._TKC), lenData)
113data
114(+1.24911904, +0.759346604), (+1.29199243, +0.813873291), (+1.23111391, +0.784803092), (+1.44666255, +0.224591851), (+1.66711354, +0.280490339), (+1.10944080, +0.357776463), (+1.64417553, +0.330640495), (+1.90659308, +0.858695507E-1), (+1.18577337, +0.754107296), (+1.55002618, +0.842631578), (+1.41862249, +0.765983939), (+1.65896153, +0.966683567), (+1.83392739, +0.954725862), (+1.60836327, +0.295341015E-1), (+1.90368962, +0.890314877), (+1.61919343, +0.430764496), (+1.90281069, +0.862413764), (+1.70729649, +0.898456454), (+1.33360839, +0.673749804), (+1.20840645, +0.588285446), (+1.21096683, +0.627702892), (+1.28788686, +0.605283558), (+1.34914160, +0.491455972), (+1.13924146, +0.317627192), (+1.30492246, +0.278793693), (+1.89244926, +0.596034348), (+1.83241630, +0.238792300), (+1.06019449, +0.575090289), (+1.66856515, +0.996734202), (+1.88007176, +0.296393454), (+1.01433086, +0.995915055), (+1.35952449, +0.225192606), (+1.56428742, +0.109285653), (+1.26882803, +0.476378381), (+1.10666299, +0.873525977), (+1.31214404, +0.831088424), (+1.70859170, +0.725056529), (+1.88534570, +0.569926262), (+1.50771940, +0.964242935), (+1.90027535, +0.216925442), (+1.17709386, +0.328605413), (+1.74551499, +0.229756296), (+1.36819911, +0.155273616), (+1.70671153, +0.328200340), (+1.73015821, +0.719621658), (+1.79838443, +0.449505508), (+1.08871627, +0.378537714), (+1.37322497, +0.680306494), (+1.97136438, +0.177314699), (+1.98874617, +0.735731542), (+1.86918974, +0.481148362), (+1.60759211, +0.137621224), (+1.81256139, +0.228771448), (+1.61021447, +0.589978695), (+1.91944671, +0.614286184), (+1.09313798, +0.576257944), (+1.28915405, +0.348247588), (+1.58724177, +0.970139205), (+1.64637899, +0.407550812), (+1.38223934, +0.966984034E-2), (+1.36540914, +0.388863087), (+1.96090174, +0.414497852E-1), (+1.38981616, +0.902015448), (+1.59574330, +0.881099701E-1), (+1.38760638, +0.583599806E-1), (+1.34131980, +0.742899477), (+1.51056230, +0.373663843), (+1.67387354, +0.211729884), (+1.79125977, +0.999856591), (+1.47659588, +0.336592197E-1), (+1.08805418, +0.192368627E-1), (+1.25629473, +0.353291392), (+1.74426770, +0.578628302), (+1.41036201, +0.169356346), (+1.43750310, +0.573294222), (+1.90894043, +0.887255430), (+1.50566852, +0.520831347E-1), (+1.54211402, +0.150596321), (+1.51236391, +0.392755806), (+1.72398639, +0.593448281E-1), (+1.16414547, +0.542320490), (+1.10033715, +0.799876451), (+1.66418600, +0.615896761), (+1.88375545, +0.987690091), (+1.50189328, +0.311064541), (+1.88035798, +0.273822844), (+1.64981103, +0.386064470), (+1.77262568, +0.988786340), (+1.56926966, +0.983290613), (+1.74002457, +0.594012499), (+1.53900218, +0.563903630), (+1.35560513, +0.738333285), (+1.01444983, +0.967131853)
115
116atad = getFFTR(getFFTF(data))
117atad = atad / size(atad)
118atad
119(+1.24911904, +0.759346604), (+1.29199231, +0.813873231), (+1.23111379, +0.784803152), (+1.44666243, +0.224591896), (+1.66711354, +0.280490220), (+1.10944068, +0.357776284), (+1.64417529, +0.330640286), (+1.90659308, +0.858693272E-1), (+1.18577349, +0.754107296), (+1.55002630, +0.842631578), (+1.41862249, +0.765983999), (+1.65896165, +0.966683567), (+1.83392763, +0.954725862), (+1.60836351, +0.295340270E-1), (+1.90368986, +0.890314877), (+1.61919367, +0.430764467), (+1.90281057, +0.862413883), (+1.70729637, +0.898456514), (+1.33360863, +0.673749804), (+1.20840633, +0.588285565), (+1.21096683, +0.627702951), (+1.28788698, +0.605283499), (+1.34914172, +0.491455972), (+1.13924158, +0.317627251), (+1.30492258, +0.278793752), (+1.89244926, +0.596034467), (+1.83241630, +0.238792419), (+1.06019449, +0.575090170), (+1.66856492, +0.996734202), (+1.88007176, +0.296393603), (+1.01433098, +0.995914936), (+1.35952449, +0.225192621), (+1.56428719, +0.109285533), (+1.26882815, +0.476378322), (+1.10666323, +0.873525739), (+1.31214404, +0.831088185), (+1.70859170, +0.725056350), (+1.88534570, +0.569925964), (+1.50771952, +0.964242756), (+1.90027547, +0.216925129), (+1.17709398, +0.328605294), (+1.74551511, +0.229756251), (+1.36819923, +0.155273527), (+1.70671165, +0.328200161), (+1.73015821, +0.719621480), (+1.79838443, +0.449505538), (+1.08871627, +0.378537685), (+1.37322521, +0.680306315), (+1.97136426, +0.177314609), (+1.98874617, +0.735731363), (+1.86918974, +0.481148243), (+1.60759187, +0.137621224), (+1.81256139, +0.228771299), (+1.61021471, +0.589978576), (+1.91944647, +0.614285946), (+1.09313822, +0.576257825), (+1.28915405, +0.348247468), (+1.58724165, +0.970139027), (+1.64637899, +0.407550752), (+1.38223958, +0.966978073E-2), (+1.36540902, +0.388862908), (+1.96090174, +0.414497852E-1), (+1.38981605, +0.902015388), (+1.59574318, +0.881097689E-1), (+1.38760638, +0.583600104E-1), (+1.34131992, +0.742899358), (+1.51056254, +0.373663813), (+1.67387342, +0.211729690), (+1.79126000, +0.999856412), (+1.47659600, +0.336592942E-1), (+1.08805394, +0.192368627E-1), (+1.25629497, +0.353291333), (+1.74426782, +0.578628242), (+1.41036212, +0.169356346), (+1.43750286, +0.573294222), (+1.90894043, +0.887255490), (+1.50566864, +0.520832241E-1), (+1.54211390, +0.150596455), (+1.51236367, +0.392755926), (+1.72398627, +0.593449473E-1), (+1.16414547, +0.542320490), (+1.10033727, +0.799876273), (+1.66418600, +0.615896583), (+1.88375533, +0.987689853), (+1.50189352, +0.311064422), (+1.88035810, +0.273822963), (+1.64981091, +0.386064410), (+1.77262580, +0.988786161), (+1.56926966, +0.983290613), (+1.74002445, +0.594012558), (+1.53900218, +0.563903570), (+1.35560513, +0.738333344), (+1.01444972, +0.967131913), (+0.00000000, -0.119209290E-6), (+0.238418579E-6, +0.119209290E-6), (+0.357627869E-6, +0.745058060E-7), (+0.596046448E-7, +0.182539225E-6), (-0.596046448E-7, +0.149011612E-7), (+0.178813934E-6, +0.238418579E-6), (+0.178813934E-6, +0.327825546E-6), (+0.298023224E-6, +0.536441803E-6), (+0.596046448E-7, +0.357627869E-6), (+0.00000000, +0.417232513E-6), (+0.417232513E-6, +0.402331352E-6), (+0.00000000, +0.745058060E-7), (+0.00000000, +0.149011612E-7), (+0.00000000, +0.894069672E-7), (+0.238418579E-6, +0.119209290E-6), (+0.119209290E-6, +0.178813934E-6), (+0.596046448E-7, +0.298023224E-7), (+0.596046448E-7, +0.268220901E-6), (+0.178813934E-6, +0.238418579E-6), (+0.238418579E-6, +0.894069672E-7), (+0.298023224E-6, +0.894069672E-7), (+0.417232513E-6, +0.298023224E-6), (+0.178813934E-6, +0.119209290E-6), (+0.357627869E-6, +0.208616257E-6), (+0.298023224E-6, +0.149011612E-6), (+0.238418579E-6, +0.208616257E-6), (+0.238418579E-6, +0.149011612E-6), (+0.476837158E-6, -0.104308128E-6), (+0.298023224E-6, -0.268220901E-6), (+0.357627869E-6, +0.745058060E-7), (+0.536441803E-6, -0.119209290E-6), (+0.119209290E-6, +0.447034836E-7), (+0.178813934E-6, -0.894069672E-7), (-0.596046448E-7, +0.178813934E-6), (+0.298023224E-6, +0.163912773E-6)
120reltol = sqrt(epsilon(1._TKC))
121reltol
122+0.345266977E-3
123isClose(data, atad(1:size(data)), reltol = reltol)
124T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T
125call setAsserted(all(isClose(data, atad(1:size(data)), reltol = reltol)))
126
127lenData = getUnifRand(5, 100)
128lenData
129+97
130data = 1._TKC + getUnifRand((0._TKC, 0._TKC), (1._TKC, 1._TKC), lenData)
131data
132(+1.1847152949006845, +0.21467341798677586), (+1.3578290563870725, +0.65320528080739837), (+1.0653746669982100, +0.82078413475330414), (+1.2000873775205725, +0.24010615531987367), (+1.5834931686033131, +0.12003890106243487), (+1.1275609590987905, +0.93034950146062834), (+1.3340936665560137, +0.83644129558827329), (+1.6296254427634183, +0.82561350183317517), (+1.3881222786082485, +0.35499486583887285), (+1.3254749910972761, +0.32830931839513122), (+1.8023347054576921, +0.32748004113751561E-1), (+1.5787269911990487, +0.95909571525265569), (+1.0783334685123003, +0.74114719175628951E-1), (+1.6123223936759874, +0.45473010392504176), (+1.5921085375129074, +0.92140339532733195), (+1.2478321825966441, +0.49010026185324573E-1), (+1.7058042248253429, +0.69115995651066153), (+1.5807527534027093, +0.90314965756842802), (+1.8572716899285631, +0.36645462432242426), (+1.5317194758630137, +0.68124332152849709), (+1.6916479968750264, +0.73104575082549683), (+1.9828546455893101, +0.53462785338647645), (+1.4408359621659290, +0.26773564450838550), (+1.0667151784210049, +0.83167772492240288), (+1.7418737916931408, +0.34581315422494385), (+1.9772511737767045, +0.27157705695837586), (+1.6639773949419752, +0.49696271019352412), (+1.8157845300266502, +0.91067487010531301), (+1.2980375369847330, +0.88741053595670949), (+1.2283059591887193, +0.66833509408278524), (+1.2540385258617863, +0.21922837822073971), (+1.5994570996040074, +0.89129790232653572), (+1.9324970959562640, +0.85200050710246200), (+1.2407903058056351, +0.92051149664891518), (+1.6202085447063486, +0.18847500864758182), (+1.5127320263749717, +0.81269674250698098), (+1.4325519068500658, +0.48311223070433795), (+1.5570683501885529, +0.96167922961860886), (+1.3599262447287022, +0.32924519448085743), (+1.7656532957106894, +0.18988292988655719), (+1.4468928288989351, +0.37667501825485428), (+1.4856852322613632, +0.22627444136368713), (+1.2376873144784375, +0.85695040086908902), (+1.2766059123105191, +0.79567671554960617), (+1.3180814209024376, +0.30702319718442628), (+1.6954432065412779, +0.23467285743838584), (+1.5768755249758057, +0.44132886679221672), (+1.0222675244365202, +0.23451760729691307), (+1.0625824323916051, +0.23340263886577428), (+1.1402439576537162, +0.84361641791515163), (+1.0155872982065821, +0.14953995596053382), (+1.0087073809938056, +0.71829580506512247), (+1.4039900602292565, +0.83675722988488854), (+1.0316950243725738, +0.64168844552648263), (+1.8168999542339606, +0.68657151162472374), (+1.2000696030249638, +0.47212258383869443), (+1.0778180361292504, +0.46637968321513468), (+1.7845355546110859, +0.49870039810279998), (+1.9712484948273321, +0.83839415107089987), (+1.3691737171100180, +0.42408357843003930), (+1.8533611587425525, +0.77330597657261368), (+1.5438256605230554, +0.15644116279862319), (+1.1779475296010160, +0.85408510466907062), (+1.2134267888291215, +0.39680726747124995), (+1.6781024781316947, +0.98428503914451670), (+1.9854390163008846, +0.21259958079949604), (+1.3409317703894321, +0.35343972883838493), (+1.6867165224293510, +0.47308412899230357), (+1.9138041437791884, +0.18150672930493461), (+1.0220747315283638, +0.77747537491713103), (+1.9676881673678619, +0.33891859667163871), (+1.6446125808966743, +0.68329734331040970), (+1.0681772867466932, +0.54304459336780064), (+1.3227157993603877, +0.73993113531722277), (+1.7802648221907975, +0.22294262044903757), (+1.3093067486212360, +0.65571397251131625), (+1.5799491371157475, +0.18451458186132841), (+1.2253306650738236, +0.40473999834633279), (+1.3757376581349912, +0.68494221930070187), (+1.5422710364481127, +0.70650744671048971), (+1.7642389810559203, +0.76072071062036883), (+1.3566913185396541, +0.91420575389391168), (+1.9593481749362089, +0.79143373353215118), (+1.1579082748325276, +0.98796077057979881), (+1.2742316442849155, +0.44122987496230470), (+1.6046445566372469, +0.14303475041174674), (+1.2986763816614793, +0.88899849175614554), (+1.5930967319817309, +0.21199480012400573E-1), (+1.5346779249698261, +0.14499793652489368), (+1.6635750243853946, +0.67684653923858307), (+1.8849075143231731, +0.34448900067072830), (+1.4393390461975168, +0.29020815835505387), (+1.8339858585422988, +0.19706590246911793E-1), (+1.4009862106508952, +0.91216082046030156E-1), (+1.0487983906710423, +0.41811534811770246), (+1.2294216935240980, +0.61751223739931560), (+1.9818888126160332, +0.80473002773958513)
133
134atad = getFFTR(getFFTF(data))
135atad = atad / size(atad)
136atad
137(+1.1847152949006849, +0.21467341798677581), (+1.3578290563870730, +0.65320528080739859), (+1.0653746669982100, +0.82078413475330403), (+1.2000873775205725, +0.24010615531987417), (+1.5834931686033131, +0.12003890106243482), (+1.1275609590987909, +0.93034950146062867), (+1.3340936665560141, +0.83644129558827318), (+1.6296254427634185, +0.82561350183317528), (+1.3881222786082485, +0.35499486583887280), (+1.3254749910972767, +0.32830931839513122), (+1.8023347054576924, +0.32748004113751616E-1), (+1.5787269911990487, +0.95909571525265580), (+1.0783334685123005, +0.74114719175629160E-1), (+1.6123223936759874, +0.45473010392504198), (+1.5921085375129074, +0.92140339532733206), (+1.2478321825966443, +0.49010026185324795E-1), (+1.7058042248253431, +0.69115995651066164), (+1.5807527534027095, +0.90314965756842802), (+1.8572716899285631, +0.36645462432242432), (+1.5317194758630135, +0.68124332152849698), (+1.6916479968750262, +0.73104575082549661), (+1.9828546455893103, +0.53462785338647634), (+1.4408359621659297, +0.26773564450838533), (+1.0667151784210049, +0.83167772492240277), (+1.7418737916931408, +0.34581315422494385), (+1.9772511737767049, +0.27157705695837597), (+1.6639773949419749, +0.49696271019352445), (+1.8157845300266502, +0.91067487010531267), (+1.2980375369847330, +0.88741053595670927), (+1.2283059591887195, +0.66833509408278502), (+1.2540385258617868, +0.21922837822073982), (+1.5994570996040067, +0.89129790232653550), (+1.9324970959562642, +0.85200050710246211), (+1.2407903058056355, +0.92051149664891507), (+1.6202085447063488, +0.18847500864758171), (+1.5127320263749717, +0.81269674250698043), (+1.4325519068500658, +0.48311223070433779), (+1.5570683501885534, +0.96167922961860863), (+1.3599262447287024, +0.32924519448085721), (+1.7656532957106896, +0.18988292988655664), (+1.4468928288989353, +0.37667501825485405), (+1.4856852322613636, +0.22627444136368710), (+1.2376873144784377, +0.85695040086908914), (+1.2766059123105191, +0.79567671554960573), (+1.3180814209024376, +0.30702319718442594), (+1.6954432065412783, +0.23467285743838562), (+1.5768755249758060, +0.44132886679221667), (+1.0222675244365202, +0.23451760729691301), (+1.0625824323916055, +0.23340263886577436), (+1.1402439576537164, +0.84361641791515152), (+1.0155872982065826, +0.14953995596053382), (+1.0087073809938056, +0.71829580506512203), (+1.4039900602292570, +0.83675722988488799), (+1.0316950243725742, +0.64168844552648241), (+1.8168999542339610, +0.68657151162472330), (+1.2000696030249638, +0.47212258383869399), (+1.0778180361292504, +0.46637968321513446), (+1.7845355546110861, +0.49870039810280009), (+1.9712484948273323, +0.83839415107089954), (+1.3691737171100180, +0.42408357843003902), (+1.8533611587425525, +0.77330597657261313), (+1.5438256605230554, +0.15644116279862308), (+1.1779475296010162, +0.85408510466907028), (+1.2134267888291217, +0.39680726747124989), (+1.6781024781316951, +0.98428503914451682), (+1.9854390163008850, +0.21259958079949623), (+1.3409317703894323, +0.35343972883838520), (+1.6867165224293510, +0.47308412899230384), (+1.9138041437791886, +0.18150672930493450), (+1.0220747315283645, +0.77747537491713081), (+1.9676881673678621, +0.33891859667163876), (+1.6446125808966745, +0.68329734331040959), (+1.0681772867466934, +0.54304459336780075), (+1.3227157993603880, +0.73993113531722299), (+1.7802648221907975, +0.22294262044903745), (+1.3093067486212360, +0.65571397251131613), (+1.5799491371157475, +0.18451458186132835), (+1.2253306650738238, +0.40473999834633279), (+1.3757376581349912, +0.68494221930070187), (+1.5422710364481125, +0.70650744671048959), (+1.7642389810559205, +0.76072071062036861), (+1.3566913185396539, +0.91420575389391179), (+1.9593481749362089, +0.79143373353215130), (+1.1579082748325276, +0.98796077057979859), (+1.2742316442849160, +0.44122987496230448), (+1.6046445566372471, +0.14303475041174671), (+1.2986763816614799, +0.88899849175614509), (+1.5930967319817309, +0.21199480012400407E-1), (+1.5346779249698259, +0.14499793652489357), (+1.6635750243853944, +0.67684653923858318), (+1.8849075143231724, +0.34448900067072830), (+1.4393390461975171, +0.29020815835505376), (+1.8339858585422988, +0.19706590246911626E-1), (+1.4009862106508957, +0.91216082046029878E-1), (+1.0487983906710425, +0.41811534811770223), (+1.2294216935240978, +0.61751223739931493), (+1.9818888126160334, +0.80473002773958502), (+0.11102230246251565E-15, -0.22204460492503131E-15), (+0.11102230246251565E-15, +0.11102230246251565E-15), (-0.11102230246251565E-15, -0.44408920985006262E-15), (+0.33306690738754696E-15, +0.22204460492503131E-15), (+0.11102230246251565E-15, -0.11102230246251565E-15), (+0.0000000000000000, +0.33306690738754696E-15), (+0.33306690738754696E-15, -0.11102230246251565E-15), (+0.0000000000000000, -0.11102230246251565E-15), (+0.0000000000000000, -0.36082248300317588E-15), (+0.0000000000000000, -0.16653345369377348E-15), (-0.11102230246251565E-15, -0.33306690738754696E-15), (-0.11102230246251565E-15, +0.13877787807814457E-15), (+0.55511151231257827E-15, +0.11102230246251565E-15), (+0.22204460492503131E-15, -0.55511151231257827E-16), (+0.44408920985006262E-15, -0.55511151231257827E-16), (+0.55511151231257827E-15, +0.24980018054066022E-15), (+0.0000000000000000, +0.55511151231257827E-16), (+0.22204460492503131E-15, +0.0000000000000000), (+0.22204460492503131E-15, +0.16653345369377348E-15), (+0.44408920985006262E-15, +0.27755575615628914E-15), (+0.77715611723760958E-15, +0.33306690738754696E-15), (+0.55511151231257827E-15, +0.49960036108132044E-15), (+0.66613381477509392E-15, +0.33306690738754696E-15), (+0.11102230246251565E-15, -0.27755575615628914E-16), (+0.22204460492503131E-15, -0.27755575615628914E-16), (+0.22204460492503131E-15, +0.11102230246251565E-15), (+0.55511151231257827E-15, +0.22204460492503131E-15), (+0.33306690738754696E-15, +0.38857805861880479E-15), (+0.77715611723760958E-15, +0.44408920985006262E-15), (+0.22204460492503131E-15, +0.38857805861880479E-15), (+0.99920072216264089E-15, +0.61062266354383610E-15)
138reltol = sqrt(epsilon(1._TKC))
139reltol
140+0.14901161193847656E-7
141isClose(data, atad(1:size(data)), reltol = reltol)
142T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T
143call setAsserted(all(isClose(data, atad(1:size(data)), reltol = reltol)))
144
145lenData = getUnifRand(5, 100)
146lenData
147+49
148data = 1._TKC + getUnifRand((0._TKC, 0._TKC), (1._TKC, 1._TKC), lenData)
149data
150(+1.57895082206820090878334940069858297, +0.308482569299724791312982041385913557), (+1.30544722772596868364063267565551659, +0.281957433943501154261889166878117212), (+1.97114232012557694495575434367952300, +0.222392006432771716396747635505585413), (+1.86699378861410385570703166523014583, +0.465211855884608968006598569313827381), (+1.51318321755707591222194074634679517, +0.580123126136185582431998811154321508), (+1.54193243766038454213137283835453564, +0.959170537253446058774125185724076442), (+1.79888113463806625286726975647277340, +0.816348892957726021596613195527074843), (+1.89312652383294487775226612070430521, +0.793834050436423952113010607825654389), (+1.33722508080128543199621348387139129, +0.366576297815534267511372227864318058), (+1.25815576949991331556895107359967769, +0.783209476446135198056711370966617869), (+1.71546545888737889555079851104981807, +0.836886207774738824230485660992476505), (+1.47650227729374227262203657587919220, +0.900325013120853408688384850330543069E-1), (+1.85635261182363125752637368151910976, +0.922109998083082698015410528990097382), (+1.86728236300314917618615051958016373, +0.408563655131725709611487425653815179), (+1.40304900571009573083234479057562492, +0.723879382087011361302085397457374313), (+1.68137099260274974747678537662926821, +0.921078661247353561578493184966320968E-2), (+1.47247108110913557056715816700812631, +0.972552622671002233831084298964563561), (+1.61851095833807817010703007423711951, +0.686183494839194976019058563157619137E-2), (+1.03014359299830774969765504512050755, +0.721020683165552394258154475128597898), (+1.66363094239158875482213885498627072, +0.949041879155550172915933227092113765), (+1.08627523190232801416008419615559589, +0.905710950985357644025733224813511849), (+1.54273846647062131304266236672952775, +0.562388006808957775827988783808964678), (+1.32073300630501557530905874121876124, +0.408751780557694769353055711445463610), (+1.23066903664151430742114012511492706, +0.358279149777114055882707712661971899), (+1.94772874819233047317643725669020086, +0.713551499016951597629692254208763791E-1), (+1.46041013964139602885239982506512536, +0.359154608309743830130304721382518208), (+1.55005376738202845311855016277297120, +0.388984876335107658897148602487937876), (+1.15449739644220470419144636436451214, +0.729934566332697942477242604326081747), (+1.84303126626867503243574373725027196, +0.254251023747785039247924618441579511), (+1.79619707564707407048031027327099928, +0.907476112827359860664962195527525862), (+1.33006387395664275309318830450533208, +0.812023853860296972880159776281492435), (+1.97371550527863651841003157060266927, +0.947313017018322097128479771274378347), (+1.87419487041101973876763424759031874, +0.724887394235918355337483190223358559), (+1.65817621574876455196223431519820220, +0.777680777281089660630969747236425898), (+1.72675228632295307127809142038983948, +0.109078307695302327829926314594412118), (+1.25474579012838865855889982171742704, +0.527951700468247410732957755150266598), (+1.62976205536831291539352791075749042, +0.356582013596378347601213315504640878E-1), (+1.74921594128973959199969117411956256, +0.323729840019365203876257158695953610), (+1.08319940283914025267350617025116307, +0.336452059700707099816829627071063946), (+1.65325401983089316464073733246934229, +0.668087657781189489359149838201818953), (+1.15151477236371278966290881008199849, +0.550097364003891025302169987504793545), (+1.43629452464512923272822752405787915, +0.335673519486888243964263813764142815), (+1.93835532866755467802715268119707953, +0.309362180083652498435782338645742909), (+1.36006798490583990323415688413840539, +0.600900932337933826396957838002696124), (+1.74281414267485920057080777042554202, +0.441136775719679304533974529634158316), (+1.42924299841949124668817466754451178, +0.329214055213337239726047003945381313), (+1.99845534244049664723585995284380115, +0.458141991841450626999633699551431045), (+1.49249602084457984078553902309624629, +0.703400640895034881911610690212497130), (+1.90774876925961266207475409923782140, +0.972042644191945922473285461091995342)
151
152atad = getFFTR(getFFTF(data))
153atad = atad / size(atad)
154atad
155(+1.57895082206820090878334940069858335, +0.308482569299724791312982041385913461), (+1.30544722772596868364063267565551659, +0.281957433943501154261889166878117212), (+1.97114232012557694495575434367952280, +0.222392006432771716396747635505585558), (+1.86699378861410385570703166523014583, +0.465211855884608968006598569313827381), (+1.51318321755707591222194074634679555, +0.580123126136185582431998811154321604), (+1.54193243766038454213137283835453583, +0.959170537253446058774125185724076346), (+1.79888113463806625286726975647277359, +0.816348892957726021596613195527074746), (+1.89312652383294487775226612070430540, +0.793834050436423952113010607825654196), (+1.33722508080128543199621348387139109, +0.366576297815534267511372227864318154), (+1.25815576949991331556895107359967788, +0.783209476446135198056711370966617966), (+1.71546545888737889555079851104981769, +0.836886207774738824230485660992476505), (+1.47650227729374227262203657587919240, +0.900325013120853408688384850330547884E-1), (+1.85635261182363125752637368151911014, +0.922109998083082698015410528990097767), (+1.86728236300314917618615051958016411, +0.408563655131725709611487425653815468), (+1.40304900571009573083234479057562473, +0.723879382087011361302085397457374409), (+1.68137099260274974747678537662926783, +0.921078661247353561578493184966369117E-2), (+1.47247108110913557056715816700812670, +0.972552622671002233831084298964563657), (+1.61851095833807817010703007423711951, +0.686183494839194976019058563157627563E-2), (+1.03014359299830774969765504512050774, +0.721020683165552394258154475128597898), (+1.66363094239158875482213885498627072, +0.949041879155550172915933227092113862), (+1.08627523190232801416008419615559570, +0.905710950985357644025733224813511945), (+1.54273846647062131304266236672952813, +0.562388006808957775827988783808964678), (+1.32073300630501557530905874121876124, +0.408751780557694769353055711445463562), (+1.23066903664151430742114012511492745, +0.358279149777114055882707712661971610), (+1.94772874819233047317643725669020086, +0.713551499016951597629692254208764272E-1), (+1.46041013964139602885239982506512555, +0.359154608309743830130304721382518352), (+1.55005376738202845311855016277297101, +0.388984876335107658897148602487937876), (+1.15449739644220470419144636436451214, +0.729934566332697942477242604326081747), (+1.84303126626867503243574373725027196, +0.254251023747785039247924618441579319), (+1.79619707564707407048031027327099967, +0.907476112827359860664962195527525862), (+1.33006387395664275309318830450533189, +0.812023853860296972880159776281492339), (+1.97371550527863651841003157060266927, +0.947313017018322097128479771274378251), (+1.87419487041101973876763424759031893, +0.724887394235918355337483190223358366), (+1.65817621574876455196223431519820220, +0.777680777281089660630969747236426090), (+1.72675228632295307127809142038983928, +0.109078307695302327829926314594412215), (+1.25474579012838865855889982171742742, +0.527951700468247410732957755150266694), (+1.62976205536831291539352791075749062, +0.356582013596378347601213315504640397E-1), (+1.74921594128973959199969117411956275, +0.323729840019365203876257158695953802), (+1.08319940283914025267350617025116307, +0.336452059700707099816829627071063753), (+1.65325401983089316464073733246934248, +0.668087657781189489359149838201819146), (+1.15151477236371278966290881008199849, +0.550097364003891025302169987504793448), (+1.43629452464512923272822752405787915, +0.335673519486888243964263813764142671), (+1.93835532866755467802715268119707915, +0.309362180083652498435782338645742909), (+1.36006798490583990323415688413840539, +0.600900932337933826396957838002696028), (+1.74281414267485920057080777042554202, +0.441136775719679304533974529634158653), (+1.42924299841949124668817466754451178, +0.329214055213337239726047003945381602), (+1.99845534244049664723585995284380134, +0.458141991841450626999633699551431141), (+1.49249602084457984078553902309624629, +0.703400640895034881911610690212497034), (+1.90774876925961266207475409923782160, +0.972042644191945922473285461091995438), (-0.962964972193617926527988971292463659E-34, +0.120370621524202240815998621411557957E-34), (+0.385185988877447170611195588516985464E-33, +0.00000000000000000000000000000000000), (+0.288889491658085377958396691387739098E-33, +0.481482486096808963263994485646231830E-34), (-0.577778983316170755916793382775478196E-33, +0.962964972193617926527988971292463659E-34), (-0.481482486096808963263994485646231830E-33, +0.00000000000000000000000000000000000), (-0.192592994438723585305597794258492732E-33, +0.288889491658085377958396691387739098E-33), (-0.192592994438723585305597794258492732E-33, +0.168518870133883137142398069976181140E-33), (+0.385185988877447170611195588516985464E-33, +0.481482486096808963263994485646231830E-34), (-0.288889491658085377958396691387739098E-33, -0.144444745829042688979198345693869549E-33), (+0.00000000000000000000000000000000000, +0.385185988877447170611195588516985464E-33), (+0.481482486096808963263994485646231830E-33, +0.962964972193617926527988971292463659E-34), (+0.962964972193617926527988971292463659E-34, -0.192592994438723585305597794258492732E-33), (+0.962964972193617926527988971292463659E-34, -0.481482486096808963263994485646231830E-33), (-0.192592994438723585305597794258492732E-33, -0.240741243048404481631997242823115915E-33), (+0.962964972193617926527988971292463659E-34, -0.962964972193617926527988971292463659E-34)
156reltol = sqrt(epsilon(1._TKC))
157reltol
158+0.138777878078144567552953958511352539E-16
159isClose(data, atad(1:size(data)), reltol = reltol)
160T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T
161call setAsserted(all(isClose(data, atad(1:size(data)), reltol = reltol)))
162
163lenData = getUnifRand(5, 100)
164lenData
165+91
166data = 1._TKC + getUnifRand(0._TKC, 1._TKC, lenData)
167data
168+1.38091278, +1.26670265, +1.63673770, +1.78869486, +1.16288674, +1.67587185, +1.64385903, +1.79863095, +1.70363712, +1.02455807, +1.84156895, +1.17303693, +1.86677980, +1.32045305, +1.42923069, +1.97866416, +1.17290199, +1.12806153, +1.66337776, +1.11898839, +1.49416971, +1.30233920, +1.20278525, +1.72135186, +1.42290020, +1.11793017, +1.59715557, +1.53116417, +1.05908108, +1.65702713, +1.17734671, +1.01181090, +1.33277202, +1.94274044, +1.57754397, +1.86223316, +1.18917787, +1.46506250, +1.36277604, +1.69898248, +1.95034409, +1.99565828, +1.70911956, +1.39411235, +1.25070953, +1.00442028, +1.84288359, +1.26464808, +1.68103862, +1.07509971, +1.80108881, +1.39493680, +1.03924346, +1.18684697, +1.71520376, +1.18948543, +1.75451064, +1.57572246, +1.19699335, +1.62439775, +1.17208970, +1.68016338, +1.24697959, +1.23392296, +1.34695268, +1.95334530, +1.37544632, +1.83540761, +1.13925540, +1.55666351, +1.59919584, +1.46424842, +1.61403513, +1.66673064, +1.46130419, +1.29266405, +1.59715199, +1.09861636, +1.56026292, +1.12560225, +1.01038122, +1.26834512, +1.34965801, +1.84768593, +1.22381926, +1.67132592, +1.16923237, +1.71804452, +1.50324380, +1.00115514, +1.61407220
169
170atad = getFFTR(getFFTF(data))
171atad = atad * 2 / size(atad)
172atad
173+1.38091278, +1.26670218, +1.63673759, +1.78869498, +1.16288662, +1.67587161, +1.64385843, +1.79863060, +1.70363712, +1.02455759, +1.84156895, +1.17303705, +1.86677980, +1.32045269, +1.42923093, +1.97866392, +1.17290211, +1.12806129, +1.66337752, +1.11898851, +1.49416935, +1.30233884, +1.20278513, +1.72135186, +1.42290044, +1.11793005, +1.59715533, +1.53116417, +1.05908108, +1.65702701, +1.17734647, +1.01181102, +1.33277225, +1.94274020, +1.57754374, +1.86223316, +1.18917799, +1.46506226, +1.36277580, +1.69898224, +1.95034409, +1.99565792, +1.70911932, +1.39411235, +1.25070977, +1.00442004, +1.84288347, +1.26464808, +1.68103838, +1.07509959, +1.80108833, +1.39493680, +1.03924334, +1.18684661, +1.71520340, +1.18948555, +1.75451052, +1.57572222, +1.19699287, +1.62439740, +1.17208934, +1.68016315, +1.24697936, +1.23392272, +1.34695268, +1.95334530, +1.37544632, +1.83540761, +1.13925552, +1.55666327, +1.59919596, +1.46424830, +1.61403513, +1.66673040, +1.46130371, +1.29266429, +1.59715199, +1.09861612, +1.56026292, +1.12560225, +1.01038098, +1.26834488, +1.34965777, +1.84768605, +1.22381914, +1.67132568, +1.16923249, +1.71804452, +1.50324368, +1.00115502, +1.61407232, +0.476837158E-6, -0.596046448E-7, +0.596046448E-7, +0.298023224E-6, +0.298023224E-6, -0.178813934E-6, +0.119209290E-6, +0.238418579E-6, +0.417232513E-6, -0.596046448E-7, +0.357627869E-6, -0.178813934E-6, +0.536441803E-6, +0.596046448E-7, -0.596046448E-7, +0.00000000, +0.476837158E-6, +0.596046448E-7, -0.178813934E-6, +0.119209290E-6, +0.238418579E-6, +0.178813934E-6, -0.119209290E-6, +0.00000000, +0.417232513E-6, +0.00000000, -0.119209290E-6, +0.119209290E-6, +0.536441803E-6, +0.357627869E-6, -0.238418579E-6, +0.178813934E-6, +0.357627869E-6, +0.119209290E-6, -0.596046448E-7, -0.119209290E-6, +0.178813934E-6
174reltol = sqrt(epsilon(1._TKC))
175reltol
176+0.345266977E-3
177isClose(data, atad(1:size(data)), reltol = reltol)
178T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T
179call setAsserted(all(isClose(data, atad(1:size(data)), reltol = reltol)))
180
181lenData = getUnifRand(5, 100)
182lenData
183+43
184data = 1._TKC + getUnifRand(0._TKC, 1._TKC, lenData)
185data
186+1.4565089448058197, +1.0766047336988076, +1.7350275319781665, +1.3725552552594897, +1.8373102710464193, +1.9525863794994698, +1.8625077345965710, +1.5949953385677076, +1.4457707790666972, +1.0773012103918151, +1.5655590511873854, +1.4936268548695113, +1.2012082042259096, +1.2138846622989452, +1.7821472178193003, +1.5869042052105353, +1.8013406541532486, +1.7915812500151493, +1.3110454822587068, +1.3126752083406812, +1.8273341973520443, +1.8139842359139373, +1.1794705789382305, +1.8087333505196148, +1.4952481296165312, +1.4399217621120544, +1.1627789293165851, +1.3089989371433970, +1.7353881990438582, +1.1479143252247646, +1.9761581139185473, +1.2145067317349074, +1.8811787968478866, +1.1025880048547969, +1.9554592584511581, +1.1142679523822063, +1.2935785648504803, +1.4748384430687076, +1.4229822142243307, +1.5872878728931272, +1.2769002079016005, +1.7694028537474376, +1.7477610524920018
187
188atad = getFFTR(getFFTF(data))
189atad = atad * 2 / size(atad)
190atad
191+1.4565089448058195, +1.0766047336988076, +1.7350275319781661, +1.3725552552594897, +1.8373102710464193, +1.9525863794994696, +1.8625077345965706, +1.5949953385677074, +1.4457707790666969, +1.0773012103918149, +1.5655590511873851, +1.4936268548695111, +1.2012082042259096, +1.2138846622989452, +1.7821472178193001, +1.5869042052105349, +1.8013406541532486, +1.7915812500151493, +1.3110454822587068, +1.3126752083406807, +1.8273341973520438, +1.8139842359139364, +1.1794705789382305, +1.8087333505196148, +1.4952481296165308, +1.4399217621120544, +1.1627789293165849, +1.3089989371433965, +1.7353881990438578, +1.1479143252247641, +1.9761581139185465, +1.2145067317349072, +1.8811787968478864, +1.1025880048547969, +1.9554592584511581, +1.1142679523822063, +1.2935785648504801, +1.4748384430687074, +1.4229822142243305, +1.5872878728931272, +1.2769002079016003, +1.7694028537474380, +1.7477610524920018, -0.33306690738754696E-15, +0.22204460492503131E-15, -0.11102230246251565E-15, +0.99920072216264089E-15, +0.0000000000000000, -0.33306690738754696E-15, +0.33306690738754696E-15, +0.44408920985006262E-15, +0.33306690738754696E-15, +0.22204460492503131E-15, -0.11102230246251565E-15, +0.22204460492503131E-15, +0.22204460492503131E-15, +0.0000000000000000, +0.11102230246251565E-15, +0.0000000000000000, +0.33306690738754696E-15, +0.0000000000000000, +0.11102230246251565E-15, +0.33306690738754696E-15, +0.66613381477509392E-15
192reltol = sqrt(epsilon(1._TKC))
193reltol
194+0.14901161193847656E-7
195isClose(data, atad(1:size(data)), reltol = reltol)
196T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T
197call setAsserted(all(isClose(data, atad(1:size(data)), reltol = reltol)))
198
199lenData = getUnifRand(5, 100)
200lenData
201+66
202data = 1._TKC + getUnifRand(0._TKC, 1._TKC, lenData)
203data
204+1.11085111540417511906131549748188240, +1.51671571870028019027584233135610314, +1.60789431016707389367763672255891849, +1.22771308608900414602552527920428549, +1.75624389701967435923795128741424805, +1.12497037568291997204908679934449603, +1.24436423541981092984547376725967269, +1.67009215156635481138168229898238689, +1.38887990873509928497894116159602247, +1.87630450367269898520598392787300848, +1.09619458398208093749151130524400586, +1.40774422452811678177658627412785037, +1.39678008637732638023882584209700062, +1.49370042614834532735624948150201224, +1.23209871916356582787921086354361728, +1.46653412253910450932433650598803702, +1.79129865892026094190626965087752901, +1.30933921682008152473178636632304311, +1.93149409026375730015392229730739207, +1.10491068233797761752777738852536788, +1.27606357649335332598747723994366419, +1.75684212462319138332243039840041217, +1.57029901977049221812954768128647669, +1.61832781076350602345342029018072540, +1.93408248456012029362377122351665551, +1.04331035360631014062155575294934668, +1.52132439591501390945975575760898322, +1.66548916337515785447613046795418548, +1.99593691879731493164047347239587029, +1.14860350379073460848399195744392482, +1.31876906446996227779575240779014441, +1.36692028815216634861725406330805341, +1.55501564920403668756694811568295044, +1.36656693931764974707335573871367513, +1.98728043532374521522202598430762491, +1.13061301666120858272112260820807134, +1.59776696634440665689845204105093150, +1.22474769632233058249460462526435253, +1.44332056623105825390465822331016962, +1.39469830727874974034643158552324800, +1.80301032618775416565707114978167030, +1.42085178837611153146290660676097180, +1.25829016870411781250478939187631420, +1.45266073205321992462687019827553466, +1.77351149692129922440525435123664616, +1.42918970723653103137352614253739651, +1.44367559050212613006288191966931082, +1.97145433585846333898066243094405125, +1.04659844745442642448324368189091413, +1.89128721712428317227278622164252487, +1.02596404050487163217890332672307936, +1.56854976112027578493727707212678579, +1.64580393708645946840161098116115415, +1.62915727551287501739631988862831708, +1.10896977026914267738170597040166450, +1.88519527891677734463735846263572891, +1.60993101819695581614495292282258028, +1.37893436615012451672829866682257959, +1.27117776791621953802333382236570592, +1.39173985422949113926976759850408788, +1.05401469242830397209100458241432583, +1.24274577551441775925569875620806204, +1.31881621649682145165381939312759242, +1.00320454172573047991403109185674539, +1.66253112343304264235018362856457917, +1.57318726998096725122443058357032992
205
206atad = getFFTR(getFFTF(data))
207atad = atad * 2 / size(atad)
208atad
209+1.11085111540417511906131549748188240, +1.51671571870028019027584233135610333, +1.60789431016707389367763672255891849, +1.22771308608900414602552527920428587, +1.75624389701967435923795128741424785, +1.12497037568291997204908679934449565, +1.24436423541981092984547376725967269, +1.67009215156635481138168229898238689, +1.38887990873509928497894116159602305, +1.87630450367269898520598392787300906, +1.09619458398208093749151130524400663, +1.40774422452811678177658627412785075, +1.39678008637732638023882584209700062, +1.49370042614834532735624948150201224, +1.23209871916356582787921086354361748, +1.46653412253910450932433650598803741, +1.79129865892026094190626965087752862, +1.30933921682008152473178636632304330, +1.93149409026375730015392229730739207, +1.10491068233797761752777738852536827, +1.27606357649335332598747723994366400, +1.75684212462319138332243039840041236, +1.57029901977049221812954768128647669, +1.61832781076350602345342029018072579, +1.93408248456012029362377122351665589, +1.04331035360631014062155575294934668, +1.52132439591501390945975575760898361, +1.66548916337515785447613046795418606, +1.99593691879731493164047347239587087, +1.14860350379073460848399195744392444, +1.31876906446996227779575240779014441, +1.36692028815216634861725406330805360, +1.55501564920403668756694811568295063, +1.36656693931764974707335573871367513, +1.98728043532374521522202598430762510, +1.13061301666120858272112260820807115, +1.59776696634440665689845204105093170, +1.22474769632233058249460462526435214, +1.44332056623105825390465822331016962, +1.39469830727874974034643158552324800, +1.80301032618775416565707114978167068, +1.42085178837611153146290660676097219, +1.25829016870411781250478939187631458, +1.45266073205321992462687019827553543, +1.77351149692129922440525435123664674, +1.42918970723653103137352614253739651, +1.44367559050212613006288191966931082, +1.97145433585846333898066243094405164, +1.04659844745442642448324368189091413, +1.89128721712428317227278622164252467, +1.02596404050487163217890332672307936, +1.56854976112027578493727707212678559, +1.64580393708645946840161098116115473, +1.62915727551287501739631988862831708, +1.10896977026914267738170597040166470, +1.88519527891677734463735846263572891, +1.60993101819695581614495292282257990, +1.37893436615012451672829866682257997, +1.27117776791621953802333382236570592, +1.39173985422949113926976759850408826, +1.05401469242830397209100458241432621, +1.24274577551441775925569875620806223, +1.31881621649682145165381939312759261, +1.00320454172573047991403109185674558, +1.66253112343304264235018362856457917, +1.57318726998096725122443058357033030, -0.962964972193617926527988971292463659E-34, -0.288889491658085377958396691387739098E-33, +0.00000000000000000000000000000000000, +0.962964972193617926527988971292463659E-34, +0.00000000000000000000000000000000000, -0.288889491658085377958396691387739098E-33, -0.481482486096808963263994485646231830E-33, +0.00000000000000000000000000000000000, -0.770371977754894341222391177033970927E-33, -0.288889491658085377958396691387739098E-33, -0.288889491658085377958396691387739098E-33, -0.288889491658085377958396691387739098E-33, -0.288889491658085377958396691387739098E-33, -0.192592994438723585305597794258492732E-33, +0.192592994438723585305597794258492732E-33, -0.385185988877447170611195588516985464E-33, -0.962964972193617926527988971292463659E-34, -0.385185988877447170611195588516985464E-33, +0.00000000000000000000000000000000000, +0.192592994438723585305597794258492732E-33, -0.962964972193617926527988971292463659E-34, -0.962964972193617926527988971292463659E-34, -0.962964972193617926527988971292463659E-34, -0.962964972193617926527988971292463659E-34, -0.288889491658085377958396691387739098E-33, -0.288889491658085377958396691387739098E-33, -0.577778983316170755916793382775478196E-33, +0.962964972193617926527988971292463659E-34, -0.962964972193617926527988971292463659E-34, +0.00000000000000000000000000000000000, -0.192592994438723585305597794258492732E-33, +0.385185988877447170611195588516985464E-33, -0.962964972193617926527988971292463659E-34, +0.192592994438723585305597794258492732E-33, -0.385185988877447170611195588516985464E-33, +0.288889491658085377958396691387739098E-33, -0.192592994438723585305597794258492732E-33, +0.00000000000000000000000000000000000, -0.962964972193617926527988971292463659E-34, -0.192592994438723585305597794258492732E-33, -0.288889491658085377958396691387739098E-33, -0.481482486096808963263994485646231830E-33, -0.577778983316170755916793382775478196E-33, -0.962964972193617926527988971292463659E-34, -0.192592994438723585305597794258492732E-33, -0.962964972193617926527988971292463659E-34, -0.962964972193617926527988971292463659E-34, +0.192592994438723585305597794258492732E-33, +0.00000000000000000000000000000000000, +0.192592994438723585305597794258492732E-33, -0.577778983316170755916793382775478196E-33, +0.962964972193617926527988971292463659E-34, +0.00000000000000000000000000000000000, +0.288889491658085377958396691387739098E-33, +0.481482486096808963263994485646231830E-33, -0.962964972193617926527988971292463659E-34, +0.288889491658085377958396691387739098E-33, -0.385185988877447170611195588516985464E-33, -0.288889491658085377958396691387739098E-33, -0.481482486096808963263994485646231830E-33, -0.288889491658085377958396691387739098E-33, -0.385185988877447170611195588516985464E-33
210reltol = sqrt(epsilon(1._TKC))
211reltol
212+0.138777878078144567552953958511352539E-16
213isClose(data, atad(1:size(data)), reltol = reltol)
214T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T
215call setAsserted(all(isClose(data, atad(1:size(data)), reltol = reltol)))
216
217lenData = getUnifRand(5, 100)
218lenData
219+7
220data = 1._TKC + getUnifRand((0._TKC, 0._TKC), (1._TKC, 1._TKC), lenData)
221data
222(+1.11822844, +0.654531062), (+1.73122239, +0.275818765), (+1.01538992, +0.730903864), (+1.69580865, +0.532580853), (+1.98573399, +0.623096645), (+1.15555108, +0.893440962), (+1.64363706, +0.389424264)
223
224atad = getFFTR(getFFTF(data))
225atad = atad / size(atad)
226atad
227(+1.11822832, +0.654531062), (+1.73122239, +0.275818974), (+1.01538992, +0.730903924), (+1.69580853, +0.532580614), (+1.98573387, +0.623096645), (+1.15555120, +0.893440843), (+1.64363694, +0.389424264), (+0.357627869E-6, +0.119209290E-6)
228reltol = sqrt(epsilon(1._TKC))
229reltol
230+0.345266977E-3
231isClose(data, atad(1:size(data)), reltol = reltol)
232T, T, T, T, T, T, T
233call setAsserted(all(isClose(data, atad(1:size(data)), reltol = reltol)))
234
235lenData = getUnifRand(5, 100)
236lenData
237+17
238data = 1._TKC + getUnifRand((0._TKC, 0._TKC), (1._TKC, 1._TKC), lenData)
239data
240(+1.6981062376411153, +0.14688174228958406E-1), (+1.5484523709811373, +0.32824525599313081), (+1.1809015761507036, +0.23705532904127791), (+1.1775593093063548, +0.52897741323042258), (+1.0644403950005283, +0.14844322441277691), (+1.6492640098851628, +0.51100792501616610), (+1.1016397980314179, +0.99306776720773804), (+1.4931701316534167, +0.74067726217094376), (+1.1101155169479595, +0.17299062773032903), (+1.5806848492826240, +0.13904389156949881), (+1.9306371620917928, +0.68530512894533080), (+1.9118287369191562, +0.54488663711917740E-1), (+1.4261582942636215, +0.18462682750953296), (+1.2685212280656413, +0.36425213839910042E-1), (+1.6308512476176249, +0.45179016356840851E-1), (+1.2180477082508485, +0.77719260501895115), (+1.3303560235364325, +0.83960290116274239)
241
242atad = getFFTR(getFFTF(data))
243atad = atad / size(atad)
244atad
245(+1.6981062376411151, +0.14688174228958462E-1), (+1.5484523709811371, +0.32824525599313098), (+1.1809015761507036, +0.23705532904127799), (+1.1775593093063548, +0.52897741323042258), (+1.0644403950005281, +0.14844322441277699), (+1.6492640098851625, +0.51100792501616610), (+1.1016397980314176, +0.99306776720773804), (+1.4931701316534163, +0.74067726217094365), (+1.1101155169479593, +0.17299062773032892), (+1.5806848492826238, +0.13904389156949870), (+1.9306371620917928, +0.68530512894533069), (+1.9118287369191560, +0.54488663711917740E-1), (+1.4261582942636215, +0.18462682750953274), (+1.2685212280656413, +0.36425213839909765E-1), (+1.6308512476176245, +0.45179016356840462E-1), (+1.2180477082508483, +0.77719260501895071), (+1.3303560235364322, +0.83960290116274239), (+0.11102230246251565E-15, -0.16653345369377348E-15), (+0.0000000000000000, -0.27755575615628914E-16), (+0.0000000000000000, +0.55511151231257827E-16), (+0.11102230246251565E-15, -0.83266726846886741E-16), (+0.22204460492503131E-15, +0.0000000000000000), (+0.22204460492503131E-15, -0.11102230246251565E-15), (+0.55511151231257827E-15, +0.55511151231257827E-16), (+0.0000000000000000, +0.41633363423443370E-16), (+0.22204460492503131E-15, +0.97144514654701197E-16), (+0.0000000000000000, -0.55511151231257827E-16), (+0.22204460492503131E-15, +0.11102230246251565E-15), (+0.33306690738754696E-15, +0.11102230246251565E-15), (+0.33306690738754696E-15, +0.38857805861880479E-15), (+0.33306690738754696E-15, +0.38857805861880479E-15), (+0.44408920985006262E-15, +0.33306690738754696E-15)
246reltol = sqrt(epsilon(1._TKC))
247reltol
248+0.14901161193847656E-7
249isClose(data, atad(1:size(data)), reltol = reltol)
250T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T
251call setAsserted(all(isClose(data, atad(1:size(data)), reltol = reltol)))
252
253lenData = getUnifRand(5, 100)
254lenData
255+15
256data = 1._TKC + getUnifRand((0._TKC, 0._TKC), (1._TKC, 1._TKC), lenData)
257data
258(+1.27839506990870126761779646498185922, +0.547225020302709191015421673949877979), (+1.17225971938159633326109332716923199, +0.219016071305788742457247231713747933), (+1.86919586375764717810145494496615910, +0.164998530513047404997100093251479586), (+1.64147371471714761716229049941328647, +0.212531745196190918963625953776959233E-1), (+1.43287019269970180227995621614702690, +0.444475026110980627343431508357543054), (+1.38439387520424353812622101836318866, +0.114757662140145614904628571202851634), (+1.90949835430134083509040591415399817, +0.898868841797070543042086754336505564), (+1.49380571389105632075792349200782623, +0.243527512925225167537057644666996411), (+1.17586012371786360594069588387792490, +0.816882230976621816368466744540138615), (+1.79933425407826468472297185939194815, +0.787606597215004022684360539850848305), (+1.77814526786786247130711511810343861, +0.665622412536402696266122079360332758), (+1.69010739224627953865863708135368607, +0.871414862661314955090934935716823197), (+1.01134503439111864827899658051162756, +0.768974479831212231308056370617062095), (+1.82433834899271217260269150461728236, +0.102355520920051226059292697713194056), (+1.98775642401957274840758615017296800, +0.957642518139523903842581586121509056)
259
260atad = getFFTR(getFFTF(data))
261atad = atad / size(atad)
262atad
263(+1.27839506990870126761779646498185942, +0.547225020302709191015421673949877979), (+1.17225971938159633326109332716923218, +0.219016071305788742457247231713747836), (+1.86919586375764717810145494496615910, +0.164998530513047404997100093251479586), (+1.64147371471714761716229049941328686, +0.212531745196190918963625953776959233E-1), (+1.43287019269970180227995621614702690, +0.444475026110980627343431508357543054), (+1.38439387520424353812622101836318885, +0.114757662140145614904628571202851622), (+1.90949835430134083509040591415399798, +0.898868841797070543042086754336505371), (+1.49380571389105632075792349200782701, +0.243527512925225167537057644666996243), (+1.17586012371786360594069588387792510, +0.816882230976621816368466744540138615), (+1.79933425407826468472297185939194834, +0.787606597215004022684360539850848593), (+1.77814526786786247130711511810343842, +0.665622412536402696266122079360332662), (+1.69010739224627953865863708135368646, +0.871414862661314955090934935716823486), (+1.01134503439111864827899658051162756, +0.768974479831212231308056370617062192), (+1.82433834899271217260269150461728255, +0.102355520920051226059292697713194212), (+1.98775642401957274840758615017296800, +0.957642518139523903842581586121508959), (-0.962964972193617926527988971292463659E-34, +0.722223729145213444895991728469347744E-34)
264reltol = sqrt(epsilon(1._TKC))
265reltol
266+0.138777878078144567552953958511352539E-16
267isClose(data, atad(1:size(data)), reltol = reltol)
268T, T, T, T, T, T, T, T, T, T, T, T, T, T, T
269call setAsserted(all(isClose(data, atad(1:size(data)), reltol = reltol)))
270
271lenData = getUnifRand(5, 100)
272lenData
273+54
274data = 1._TKC + getUnifRand(0._TKC, 1._TKC, lenData)
275data
276+1.65566111, +1.16979384, +1.64829707, +1.26386511, +1.70488870, +1.40406036, +1.72641695, +1.91135335, +1.27007604, +1.55602872, +1.39921379, +1.95771062, +1.84776282, +1.69135547, +1.66063619, +1.53991675, +1.05961466, +1.85157824, +1.14284301, +1.20080471, +1.42783165, +1.51344085, +1.46104026, +1.36868739, +1.05241323, +1.40059805, +1.42516720, +1.67613554, +1.14854848, +1.31719530, +1.97465670, +1.36911070, +1.37834358, +1.09428453, +1.12961721, +1.84418535, +1.29546404, +1.00526679, +1.33836174, +1.46081853, +1.78336716, +1.95199347, +1.06199288, +1.97592115, +1.03639746, +1.72954726, +1.11783147, +1.88619113, +1.49873900, +1.35820627, +1.20112240, +1.38081801, +1.68295431, +1.95974815
277
278atad = getFFTR(getFFTF(data))
279atad = atad * 2 / size(atad)
280atad
281+1.65566158, +1.16979373, +1.64829695, +1.26386499, +1.70488882, +1.40406036, +1.72641706, +1.91135299, +1.27007616, +1.55602860, +1.39921379, +1.95771062, +1.84776282, +1.69135559, +1.66063643, +1.53991675, +1.05961514, +1.85157824, +1.14284289, +1.20080471, +1.42783177, +1.51344097, +1.46104038, +1.36868727, +1.05241346, +1.40059781, +1.42516720, +1.67613542, +1.14854836, +1.31719542, +1.97465682, +1.36911047, +1.37834406, +1.09428442, +1.12961733, +1.84418535, +1.29546428, +1.00526690, +1.33836198, +1.46081865, +1.78336728, +1.95199323, +1.06199312, +1.97592127, +1.03639770, +1.72954738, +1.11783195, +1.88619113, +1.49873924, +1.35820627, +1.20112240, +1.38081789, +1.68295443, +1.95974815, +0.357627869E-6, +0.119209290E-6, +0.357627869E-6, +0.596046448E-7, +0.119209290E-6, +0.238418579E-6, +0.357627869E-6, +0.238418579E-6, +0.596046448E-7, +0.119209290E-6
282reltol = sqrt(epsilon(1._TKC))
283reltol
284+0.345266977E-3
285isClose(data, atad(1:size(data)), reltol = reltol)
286T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T
287call setAsserted(all(isClose(data, atad(1:size(data)), reltol = reltol)))
288
289lenData = getUnifRand(5, 100)
290lenData
291+47
292data = 1._TKC + getUnifRand(0._TKC, 1._TKC, lenData)
293data
294+1.0291980734937773, +1.6981181613681871, +1.7996166129832520, +1.7949728623974495, +1.5524092111793095, +1.9181481870076662, +1.9186294549865419, +1.3816486371931198, +1.5639665054770608, +1.3761397929628090, +1.4253132280439158, +1.3695304606967227, +1.6820585281443894, +1.8793151437443654, +1.1323471744707163, +1.5066682564299905, +1.8442054822266623, +1.2456041930204798, +1.5750363731947876, +1.5087787857049453, +1.6041435843619420, +1.7026020524250249, +1.6429228514992427, +1.0357571862090040, +1.1376892017672571, +1.8237469223986043, +1.1354227949870810, +1.9993776157472882, +1.4119313451508273, +1.1794291930720715, +1.8648568374342613, +1.9692721949518446, +1.5817633749158588, +1.9797223667340675, +1.9728078696787770, +1.6488665624691623, +1.1922558074511520, +1.6842883711941998, +1.6340179284320544, +1.4453167975263432, +1.9733945817003224, +1.7358838794079117, +1.5236019676149604, +1.9377004250342575, +1.2394344301805500, +1.2245676484559320, +1.4869422931341882
295
296atad = getFFTR(getFFTF(data))
297atad = atad * 2 / size(atad)
298atad
299+1.0291980734937773, +1.6981181613681871, +1.7996166129832518, +1.7949728623974497, +1.5524092111793093, +1.9181481870076662, +1.9186294549865415, +1.3816486371931200, +1.5639665054770604, +1.3761397929628092, +1.4253132280439156, +1.3695304606967229, +1.6820585281443889, +1.8793151437443656, +1.1323471744707163, +1.5066682564299905, +1.8442054822266623, +1.2456041930204798, +1.5750363731947872, +1.5087787857049451, +1.6041435843619418, +1.7026020524250249, +1.6429228514992429, +1.0357571862090040, +1.1376892017672571, +1.8237469223986045, +1.1354227949870808, +1.9993776157472876, +1.4119313451508271, +1.1794291930720715, +1.8648568374342611, +1.9692721949518437, +1.5817633749158588, +1.9797223667340680, +1.9728078696787767, +1.6488665624691621, +1.1922558074511518, +1.6842883711941998, +1.6340179284320548, +1.4453167975263430, +1.9733945817003220, +1.7358838794079119, +1.5236019676149601, +1.9377004250342573, +1.2394344301805498, +1.2245676484559322, +1.4869422931341885, +0.11102230246251565E-15, +0.11102230246251565E-15, +0.55511151231257827E-15, +0.0000000000000000, +0.0000000000000000, -0.22204460492503131E-15, -0.22204460492503131E-15, +0.66613381477509392E-15, +0.0000000000000000, -0.22204460492503131E-15, +0.55511151231257827E-15, +0.44408920985006262E-15, +0.66613381477509392E-15, +0.33306690738754696E-15, +0.33306690738754696E-15, +0.44408920985006262E-15, +0.88817841970012523E-15
300reltol = sqrt(epsilon(1._TKC))
301reltol
302+0.14901161193847656E-7
303isClose(data, atad(1:size(data)), reltol = reltol)
304T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T
305call setAsserted(all(isClose(data, atad(1:size(data)), reltol = reltol)))
306
307lenData = getUnifRand(5, 100)
308lenData
309+51
310data = 1._TKC + getUnifRand(0._TKC, 1._TKC, lenData)
311data
312+1.66575346949199971823403026883733315, +1.03721803084027767020726404706327484, +1.25731167782734775324910856014873178, +1.59127990664540918356792218193011394, +1.11158050320838636259101687278511579, +1.79153559129945318208958948355011960, +1.01500095709191500028298686782228655, +1.63330175932491833375807192099511353, +1.63793141579634013466880584377136297, +1.58767090733798231753971437635093378, +1.50206502989592515925233195341852910, +1.86386560353530875791440197859021374, +1.13010607598115177632678903355417029, +1.63304888484208116383705680608747317, +1.42582695984720494215609678906841657, +1.58432571352767024680573369712082502, +1.80253815623895514938491351997474322, +1.30476194463164604705452959698999536, +1.01951528921150534343820015381205366, +1.75872011559163837661025257986090033, +1.45191500980442266858002378646840801, +1.27753323262263297008124285980311685, +1.67363213663690095420726698952961419, +1.13098846175200428240194333912862851, +1.48421497580531747169135097179594023, +1.25424664507290707823632745719996866, +1.83144685981590240946161405551161097, +1.30528462308540495342300621680662296, +1.40410809380151622615362249314740970, +1.09650856074203322550341163991560360, +1.92758261607618408613075473999476582, +1.23116916949437881231891713433627473, +1.01945325900777139871230563044546033, +1.46635823448440073937892332326819667, +1.82915856018710379133316098950803786, +1.12247856335700627531802629799924977, +1.77286768926939678119143279924917213, +1.30731672841922400455754837217611311, +1.53326176099792095414705151536932838, +1.05986127495423457719380722437606588, +1.78629050450114627986692587906354263, +1.00440293712859946693823542974461616, +1.29680414456560160993315613824574714, +1.35690882630093179577780909178495140, +1.08805967146477173792645575730940421, +1.54023467537108188064473924064223445, +1.64764491295922179067448272956076737, +1.18102332777799399717735314987163152, +1.65624670568853100083763105962081288, +1.60544504074545722078701515724201066, +1.82601123267646072597804990805658392
313
314atad = getFFTR(getFFTF(data))
315atad = atad * 2 / size(atad)
316atad
317+1.66575346949199971823403026883733296, +1.03721803084027767020726404706327522, +1.25731167782734775324910856014873139, +1.59127990664540918356792218193011394, +1.11158050320838636259101687278511598, +1.79153559129945318208958948355011999, +1.01500095709191500028298686782228635, +1.63330175932491833375807192099511314, +1.63793141579634013466880584377136277, +1.58767090733798231753971437635093397, +1.50206502989592515925233195341852910, +1.86386560353530875791440197859021374, +1.13010607598115177632678903355417009, +1.63304888484208116383705680608747317, +1.42582695984720494215609678906841657, +1.58432571352767024680573369712082483, +1.80253815623895514938491351997474303, +1.30476194463164604705452959698999556, +1.01951528921150534343820015381205366, +1.75872011559163837661025257986089994, +1.45191500980442266858002378646840782, +1.27753323262263297008124285980311685, +1.67363213663690095420726698952961400, +1.13098846175200428240194333912862793, +1.48421497580531747169135097179594004, +1.25424664507290707823632745719996866, +1.83144685981590240946161405551161078, +1.30528462308540495342300621680662258, +1.40410809380151622615362249314740932, +1.09650856074203322550341163991560321, +1.92758261607618408613075473999476543, +1.23116916949437881231891713433627454, +1.01945325900777139871230563044546014, +1.46635823448440073937892332326819706, +1.82915856018710379133316098950803786, +1.12247856335700627531802629799924977, +1.77286768926939678119143279924917213, +1.30731672841922400455754837217611369, +1.53326176099792095414705151536932819, +1.05986127495423457719380722437606569, +1.78629050450114627986692587906354283, +1.00440293712859946693823542974461635, +1.29680414456560160993315613824574714, +1.35690882630093179577780909178495140, +1.08805967146477173792645575730940382, +1.54023467537108188064473924064223484, +1.64764491295922179067448272956076737, +1.18102332777799399717735314987163152, +1.65624670568853100083763105962081269, +1.60544504074545722078701515724201085, +1.82601123267646072597804990805658353, +0.962964972193617926527988971292463659E-34, -0.385185988877447170611195588516985464E-33, -0.962964972193617926527988971292463659E-34, -0.192592994438723585305597794258492732E-33, -0.674075480535532548569592279904724561E-33, -0.192592994438723585305597794258492732E-33, +0.00000000000000000000000000000000000, -0.481482486096808963263994485646231830E-33, -0.481482486096808963263994485646231830E-33, -0.288889491658085377958396691387739098E-33, -0.481482486096808963263994485646231830E-33, -0.674075480535532548569592279904724561E-33, -0.481482486096808963263994485646231830E-33
318reltol = sqrt(epsilon(1._TKC))
319reltol
320+0.138777878078144567552953958511352539E-16
321isClose(data, atad(1:size(data)), reltol = reltol)
322T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T
323call setAsserted(all(isClose(data, atad(1:size(data)), reltol = reltol)))
324
325lenData = getUnifRand(5, 100)
326lenData
327+12
328data = 1._TKC + getUnifRand((0._TKC, 0._TKC), (1._TKC, 1._TKC), lenData)
329data
330(+1.57643318, +0.990336120), (+1.52501655, +0.288927019), (+1.91373110, +0.900449753E-1), (+1.26149750, +0.868839264), (+1.91632664, +0.632703185), (+1.63471222, +0.433197677), (+1.93614388, +0.972430527), (+1.11746335, +0.281108916), (+1.79997730, +0.159388721), (+1.99071884, +0.959326088), (+1.80985129, +0.716028512), (+1.37131000, +0.531579554)
331
332atad = getFFTR(getFFTF(data))
333atad = atad / size(atad)
334atad
335(+1.57643318, +0.990336061), (+1.52501655, +0.288927138), (+1.91373098, +0.900451541E-1), (+1.26149738, +0.868839204), (+1.91632664, +0.632702947), (+1.63471246, +0.433197618), (+1.93614388, +0.972430408), (+1.11746359, +0.281108856), (+1.79997730, +0.159388781), (+1.99071884, +0.959326208), (+1.80985105, +0.716028392), (+1.37130988, +0.531579316), (-0.119209290E-6, +0.238418579E-6), (+0.00000000, +0.447034836E-7), (+0.178813934E-6, +0.178813934E-6), (+0.178813934E-6, +0.447034836E-7)
336reltol = sqrt(epsilon(1._TKC))
337reltol
338+0.345266977E-3
339isClose(data, atad(1:size(data)), reltol = reltol)
340T, T, T, T, T, T, T, T, T, T, T, T
341call setAsserted(all(isClose(data, atad(1:size(data)), reltol = reltol)))
342
343lenData = getUnifRand(5, 100)
344lenData
345+31
346data = 1._TKC + getUnifRand((0._TKC, 0._TKC), (1._TKC, 1._TKC), lenData)
347data
348(+1.6246883983756542, +0.92832952596511142), (+1.1501494447506115, +0.63307420212121790), (+1.5021588590231021, +0.79190068017224780), (+1.8988684393889323, +0.14589508724794098), (+1.4104359100846802, +0.67651513182770939), (+1.3503587682207963, +0.60284716190789733), (+1.1319741970868731, +0.49504957007639794E-1), (+1.7740439338880369, +0.49419344643546048), (+1.1871570115126229, +0.48263420155271941), (+1.7530813529781826, +0.61624574226250584), (+1.9969213638641532, +0.30569506010441694), (+1.1998455581109888, +0.31435266986609578), (+1.9495580415170035, +0.75346284029738164), (+1.6509335389251207, +0.68471734928412031), (+1.6796033098000276, +0.21760725575652995), (+1.2430469513152320, +0.38766095214204044), (+1.9910457784258229, +0.70639839658856007), (+1.3768284729390194, +0.26111827296241941), (+1.9099998508807254, +0.78859962117465143), (+1.1420433699810162, +0.72064998412490588E-1), (+1.2224938500494535, +0.28313584496540878), (+1.7705456623821156, +0.85818411295152230), (+1.8669807247915915, +0.66846627141839110), (+1.8114241450845299, +0.72757778974181342), (+1.3450659612930060, +0.84527953760544094), (+1.7618189985560599, +0.37865267909812850), (+1.5783655131026948, +0.61850216750704035), (+1.4731923526155701, +0.71009535953324987), (+1.8152388206794381, +0.51219039470013983), (+1.9194787081869424, +0.93682097574737366), (+1.4585356160872003, +0.74970161361232923)
349
350atad = getFFTR(getFFTF(data))
351atad = atad / size(atad)
352atad
353(+1.6246883983756544, +0.92832952596511131), (+1.1501494447506118, +0.63307420212121790), (+1.5021588590231023, +0.79190068017224780), (+1.8988684393889317, +0.14589508724794120), (+1.4104359100846802, +0.67651513182770939), (+1.3503587682207960, +0.60284716190789711), (+1.1319741970868733, +0.49504957007639960E-1), (+1.7740439338880365, +0.49419344643546048), (+1.1871570115126231, +0.48263420155271947), (+1.7530813529781824, +0.61624574226250561), (+1.9969213638641530, +0.30569506010441710), (+1.1998455581109884, +0.31435266986609589), (+1.9495580415170033, +0.75346284029738153), (+1.6509335389251205, +0.68471734928412020), (+1.6796033098000276, +0.21760725575653023), (+1.2430469513152316, +0.38766095214204033), (+1.9910457784258233, +0.70639839658855996), (+1.3768284729390192, +0.26111827296241930), (+1.9099998508807257, +0.78859962117465143), (+1.1420433699810157, +0.72064998412490699E-1), (+1.2224938500494535, +0.28313584496540883), (+1.7705456623821150, +0.85818411295152208), (+1.8669807247915913, +0.66846627141839110), (+1.8114241450845294, +0.72757778974181342), (+1.3450659612930063, +0.84527953760544094), (+1.7618189985560597, +0.37865267909812839), (+1.5783655131026948, +0.61850216750704035), (+1.4731923526155697, +0.71009535953324976), (+1.8152388206794381, +0.51219039470013983), (+1.9194787081869420, +0.93682097574737300), (+1.4585356160872003, +0.74970161361232912), (+0.33306690738754696E-15, +0.13877787807814457E-15)
354reltol = sqrt(epsilon(1._TKC))
355reltol
356+0.14901161193847656E-7
357isClose(data, atad(1:size(data)), reltol = reltol)
358T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T
359call setAsserted(all(isClose(data, atad(1:size(data)), reltol = reltol)))
360
361lenData = getUnifRand(5, 100)
362lenData
363+47
364data = 1._TKC + getUnifRand((0._TKC, 0._TKC), (1._TKC, 1._TKC), lenData)
365data
366(+1.84938382157359080523247154063328082, +0.647534091777212406607901729757260146), (+1.83941125756815690425962933274721897, +0.527946422941329733343012830886936017E-1), (+1.73267267309116746583659461449441597, +0.842627477412771150456831495754131291), (+1.02516730324259484620876804603852606, +0.445206535233554554054488769217694537), (+1.18021093366962994726468934883699483, +0.596774436774096835553091167618175317), (+1.47621556525223020873376870666256047, +0.182981985745376862245344206756379846), (+1.47532677530502933367312355543336952, +0.447077427812626802604554981942554057), (+1.38649838819120367664433903788788339, +0.667792717939813358241481212939745002), (+1.57760184378363768170767205855557068, +0.113307517527511888393747395550298697), (+1.64776415055395845103265070499893377, +0.491273359446976886568825483698834000), (+1.12644007917641559835910051147012207, +0.114917700078388505613592864331959273), (+1.76867870708775673549576081507175808, +0.211885629762118382414251799732938184), (+1.12737806080247841965547403180114692, +0.911475815483593405662674425278675722), (+1.79867717024210195005529880275819571, +0.290304590634212324804936912790458220), (+1.13130943579164632775087269533650399, +0.471348241280870185937978841165975485), (+1.32546964888829461689002489653256598, +0.386588271896477665135522819446281403E-1), (+1.39139182281252842376053679344142667, +0.646251433038650883838933301765986548), (+1.37317693054395722098805521985889132, +0.173417633645778145268559059315207092), (+1.63826785228908570064453861576906105, +0.141442243254844416495483368996201606), (+1.70497171534870925501118371419513873, +0.379429125970718356010733763998902745), (+1.16005720797479067328093850891796425, +0.876578554533663483397165964106621411), (+1.41376530130003669232436607234285352, +0.937242915791048561352231637445358227), (+1.41500497786586281281821081717420959, +0.505777745991044924650365334508500750), (+1.20817078998730993492321788338112959, +0.216822059814953855803034753260581828), (+1.20125727170032263171588383582064283, +0.806659670308174235849084346231270480), (+1.52738695214258719477236436279784976, +0.703114210480471493081277630235766037), (+1.96763540117032251853770590163572782, +0.654896283930733885800490769005708508), (+1.13546307503377649519559777482164718, +0.343053374227510156345033894707835838), (+1.55266178023443653551803095280061331, +0.725733819192454369052029472578067410), (+1.84511691486633354088702470211330938, +0.972429881486120752694948558706214579), (+1.72522504841744549034018495581397028, +0.178632177752791248082410744974251193), (+1.68982187332606037228461996393785093, +0.682167531442774536945319093442428579), (+1.02569803400415350804314752843832409, +0.908654598008129780476562580332027122), (+1.32741097850001427712182575211397513, +0.206731410347187656679656305303046993), (+1.16171642759812714159813097512313588, +0.304700916972541737286575742431071128), (+1.72200384681423524890926868096169280, +0.555144244725733281043740637601567038), (+1.83643897407482107995607986446326117, +0.337010293535011131423549352946250571), (+1.11990475021744193202667838093309474, +0.986211297298384991752633166996443065), (+1.95347384698244823941865352905500133, +0.573611165066146352218697945420679188E-1), (+1.92039010123401896852295651993423291, +0.646015738063764516902308625797797990E-1), (+1.34371852687488395435816773206034709, +0.400729828381877836376700662388972625), (+1.15334730678742995952467104908411781, +0.398750721194763555836988126262907537), (+1.97136070415619740356354405167891763, +0.116028288128995896078369323181830346), (+1.68486583149214274346482209582392027, +0.770355255169255198719848138725751034), (+1.13781945418751770418229205892825451, +0.502151135913384069312496409143850648E-1), (+1.98320923581582147229627484161165205, +0.167716320742087481080079245388049071), (+1.90054656595341930740549662046128130, +0.989259601000041434874572248814568020)
367
368atad = getFFTR(getFFTF(data))
369atad = atad / size(atad)
370atad
371(+1.84938382157359080523247154063328102, +0.647534091777212406607901729757260243), (+1.83941125756815690425962933274721897, +0.527946422941329733343012830886935776E-1), (+1.73267267309116746583659461449441597, +0.842627477412771150456831495754131291), (+1.02516730324259484620876804603852606, +0.445206535233554554054488769217694585), (+1.18021093366962994726468934883699483, +0.596774436774096835553091167618175317), (+1.47621556525223020873376870666256047, +0.182981985745376862245344206756379461), (+1.47532677530502933367312355543336952, +0.447077427812626802604554981942553960), (+1.38649838819120367664433903788788339, +0.667792717939813358241481212939745194), (+1.57760184378363768170767205855557087, +0.113307517527511888393747395550298793), (+1.64776415055395845103265070499893358, +0.491273359446976886568825483698833952), (+1.12644007917641559835910051147012227, +0.114917700078388505613592864331959478), (+1.76867870708775673549576081507175808, +0.211885629762118382414251799732938232), (+1.12737806080247841965547403180114692, +0.911475815483593405662674425278675818), (+1.79867717024210195005529880275819590, +0.290304590634212324804936912790458509), (+1.13130943579164632775087269533650418, +0.471348241280870185937978841165975292), (+1.32546964888829461689002489653256598, +0.386588271896477665135522819446284172E-1), (+1.39139182281252842376053679344142667, +0.646251433038650883838933301765986548), (+1.37317693054395722098805521985889112, +0.173417633645778145268559059315207188), (+1.63826785228908570064453861576906105, +0.141442243254844416495483368996201437), (+1.70497171534870925501118371419513873, +0.379429125970718356010733763998902648), (+1.16005720797479067328093850891796425, +0.876578554533663483397165964106621411), (+1.41376530130003669232436607234285313, +0.937242915791048561352231637445358420), (+1.41500497786586281281821081717420921, +0.505777745991044924650365334508500654), (+1.20817078998730993492321788338112959, +0.216822059814953855803034753260581587), (+1.20125727170032263171588383582064302, +0.806659670308174235849084346231270480), (+1.52738695214258719477236436279784976, +0.703114210480471493081277630235766230), (+1.96763540117032251853770590163572763, +0.654896283930733885800490769005708412), (+1.13546307503377649519559777482164738, +0.343053374227510156345033894707835741), (+1.55266178023443653551803095280061311, +0.725733819192454369052029472578067410), (+1.84511691486633354088702470211330918, +0.972429881486120752694948558706214772), (+1.72522504841744549034018495581397048, +0.178632177752791248082410744974250952), (+1.68982187332606037228461996393785093, +0.682167531442774536945319093442428290), (+1.02569803400415350804314752843832409, +0.908654598008129780476562580332027218), (+1.32741097850001427712182575211397513, +0.206731410347187656679656305303047209), (+1.16171642759812714159813097512313608, +0.304700916972541737286575742431071273), (+1.72200384681423524890926868096169280, +0.555144244725733281043740637601566845), (+1.83643897407482107995607986446326136, +0.337010293535011131423549352946250571), (+1.11990475021744193202667838093309435, +0.986211297298384991752633166996443065), (+1.95347384698244823941865352905500133, +0.573611165066146352218697945420677262E-1), (+1.92039010123401896852295651993423310, +0.646015738063764516902308625797797990E-1), (+1.34371852687488395435816773206034728, +0.400729828381877836376700662388972625), (+1.15334730678742995952467104908411761, +0.398750721194763555836988126262907489), (+1.97136070415619740356354405167891743, +0.116028288128995896078369323181830069), (+1.68486583149214274346482209582392066, +0.770355255169255198719848138725750552), (+1.13781945418751770418229205892825432, +0.502151135913384069312496409143851611E-1), (+1.98320923581582147229627484161165185, +0.167716320742087481080079245388049023), (+1.90054656595341930740549662046128168, +0.989259601000041434874572248814568213), (+0.00000000000000000000000000000000000, -0.361111864572606722447995864234673872E-34), (-0.962964972193617926527988971292463659E-34, -0.962964972193617926527988971292463659E-34), (+0.00000000000000000000000000000000000, -0.120370621524202240815998621411557957E-33), (+0.288889491658085377958396691387739098E-33, -0.120370621524202240815998621411557957E-33), (+0.962964972193617926527988971292463659E-34, -0.962964972193617926527988971292463659E-34), (-0.481482486096808963263994485646231830E-33, -0.192592994438723585305597794258492732E-33), (-0.674075480535532548569592279904724561E-33, -0.144444745829042688979198345693869549E-33), (-0.385185988877447170611195588516985464E-33, +0.192592994438723585305597794258492732E-33), (-0.288889491658085377958396691387739098E-33, -0.144444745829042688979198345693869549E-33), (+0.192592994438723585305597794258492732E-33, +0.00000000000000000000000000000000000), (+0.00000000000000000000000000000000000, -0.481482486096808963263994485646231830E-34), (+0.288889491658085377958396691387739098E-33, +0.144444745829042688979198345693869549E-33), (+0.577778983316170755916793382775478196E-33, -0.962964972193617926527988971292463659E-34), (-0.962964972193617926527988971292463659E-34, -0.240741243048404481631997242823115915E-33), (-0.385185988877447170611195588516985464E-33, -0.481482486096808963263994485646231830E-33), (-0.192592994438723585305597794258492732E-33, -0.156481807981462913060798207835025345E-33), (+0.00000000000000000000000000000000000, +0.962964972193617926527988971292463659E-34)
372reltol = sqrt(epsilon(1._TKC))
373reltol
374+0.138777878078144567552953958511352539E-16
375isClose(data, atad(1:size(data)), reltol = reltol)
376T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T
377call setAsserted(all(isClose(data, atad(1:size(data)), reltol = reltol)))
378
379lenData = getUnifRand(5, 100)
380lenData
381+72
382data = 1._TKC + getUnifRand(0._TKC, 1._TKC, lenData)
383data
384+1.29041147, +1.79609847, +1.05737829, +1.13525224, +1.85536492, +1.45694757, +1.06955290, +1.14406204, +1.02240348, +1.08346999, +1.79956079, +1.46049452, +1.22200263, +1.96891725, +1.94340181, +1.85811257, +1.44269919, +1.67893827, +1.11822987, +1.51502740, +1.61545658, +1.04100108, +1.43423724, +1.71532059, +1.14078498, +1.81664181, +1.82753766, +1.25474918, +1.09938002, +1.19849968, +1.70501685, +1.00769663, +1.01459718, +1.96715760, +1.45223665, +1.61774039, +1.87702274, +1.36623645, +1.03099728, +1.58792686, +1.27573180, +1.52863693, +1.38691306, +1.67848718, +1.98393106, +1.35236979, +1.53935933, +1.51672316, +1.11559510, +1.05884349, +1.96072853, +1.01130056, +1.11830556, +1.43028176, +1.62355578, +1.28529167, +1.37570679, +1.51493979, +1.11649370, +1.09698796, +1.64070678, +1.89018321, +1.32841802, +1.44221950, +1.42665672, +1.83460617, +1.93856621, +1.61596859, +1.72286606, +1.22657585, +1.86008024, +1.39717245
385
386atad = getFFTR(getFFTF(data))
387atad = atad * 2 / size(atad)
388atad
389+1.29041171, +1.79609847, +1.05737829, +1.13525212, +1.85536456, +1.45694745, +1.06955242, +1.14406180, +1.02240348, +1.08347023, +1.79956079, +1.46049452, +1.22200274, +1.96891713, +1.94340158, +1.85811234, +1.44269896, +1.67893839, +1.11822975, +1.51502752, +1.61545634, +1.04100132, +1.43423700, +1.71532059, +1.14078450, +1.81664228, +1.82753730, +1.25474942, +1.09937978, +1.19850004, +1.70501685, +1.00769699, +1.01459742, +1.96715784, +1.45223677, +1.61774039, +1.87702274, +1.36623645, +1.03099751, +1.58792686, +1.27573180, +1.52863705, +1.38691330, +1.67848730, +1.98393095, +1.35236955, +1.53935957, +1.51672316, +1.11559510, +1.05884337, +1.96072865, +1.01130056, +1.11830568, +1.43028152, +1.62355578, +1.28529143, +1.37570679, +1.51493979, +1.11649346, +1.09698796, +1.64070666, +1.89018321, +1.32841802, +1.44221973, +1.42665696, +1.83460617, +1.93856621, +1.61596859, +1.72286630, +1.22657573, +1.86008072, +1.39717245, +0.119209290E-6, -0.119209290E-6, +0.178813934E-6, +0.00000000, +0.00000000, +0.119209290E-6, +0.476837158E-6, +0.596046448E-7, +0.417232513E-6, -0.298023224E-6, +0.119209290E-6, +0.00000000, +0.536441803E-6, +0.596046448E-7, +0.596046448E-6, +0.00000000, +0.536441803E-6, -0.417232513E-6, +0.178813934E-6, -0.238418579E-6, +0.238418579E-6, -0.238418579E-6, +0.298023224E-6, +0.00000000, +0.178813934E-6, -0.178813934E-6, -0.119209290E-6, +0.596046448E-7, +0.00000000, +0.178813934E-6, -0.596046448E-7, +0.178813934E-6, +0.596046448E-7, -0.119209290E-6, +0.00000000, +0.00000000, +0.119209290E-6, +0.298023224E-6, +0.238418579E-6, +0.178813934E-6, +0.238418579E-6, -0.596046448E-7, +0.596046448E-7, +0.00000000, +0.00000000, +0.238418579E-6, +0.00000000, +0.178813934E-6, +0.238418579E-6, +0.238418579E-6, +0.178813934E-6, -0.596046448E-7, +0.119209290E-6, +0.178813934E-6, +0.596046448E-7, -0.119209290E-6
390reltol = sqrt(epsilon(1._TKC))
391reltol
392+0.345266977E-3
393isClose(data, atad(1:size(data)), reltol = reltol)
394T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T
395call setAsserted(all(isClose(data, atad(1:size(data)), reltol = reltol)))
396
397lenData = getUnifRand(5, 100)
398lenData
399+42
400data = 1._TKC + getUnifRand(0._TKC, 1._TKC, lenData)
401data
402+1.6633014772696453, +1.4872557571109839, +1.1036170521502560, +1.0152582790941602, +1.1868504705181862, +1.1506305177226346, +1.6234337071965765, +1.6868490333243731, +1.0123868941353074, +1.1073000830146689, +1.9859790624373170, +1.4439145538140745, +1.4668144209925473, +1.6260641143302845, +1.9447028289091981, +1.3794339354158232, +1.6108786282468626, +1.1122151773831481, +1.7186126473237384, +1.0796054173828191, +1.7122864470606618, +1.1381722939499741, +1.6224057963605469, +1.5648364502369378, +1.4254822369461970, +1.2742350624382963, +1.1239151951390016, +1.0985796338409433, +1.0262758120175608, +1.1619143531126475, +1.1312763232864160, +1.6577797768748042, +1.7811711518402422, +1.8757832185350742, +1.6314822883275761, +1.2880632586291307, +1.2551668490745702, +1.5206031603791430, +1.2103173948130097, +1.2679622172679401, +1.9247165921231120, +1.5564612384562593
403
404atad = getFFTR(getFFTF(data))
405atad = atad * 2 / size(atad)
406atad
407+1.6633014772696453, +1.4872557571109848, +1.1036170521502560, +1.0152582790941607, +1.1868504705181862, +1.1506305177226350, +1.6234337071965761, +1.6868490333243733, +1.0123868941353076, +1.1073000830146693, +1.9859790624373168, +1.4439145538140745, +1.4668144209925471, +1.6260641143302847, +1.9447028289091972, +1.3794339354158232, +1.6108786282468626, +1.1122151773831481, +1.7186126473237384, +1.0796054173828187, +1.7122864470606616, +1.1381722939499739, +1.6224057963605467, +1.5648364502369376, +1.4254822369461966, +1.2742350624382963, +1.1239151951390016, +1.0985796338409424, +1.0262758120175604, +1.1619143531126470, +1.1312763232864160, +1.6577797768748033, +1.7811711518402422, +1.8757832185350747, +1.6314822883275761, +1.2880632586291305, +1.2551668490745704, +1.5206031603791428, +1.2103173948130095, +1.2679622172679401, +1.9247165921231117, +1.5564612384562597, +0.44408920985006262E-15, -0.11102230246251565E-15, +0.22204460492503131E-15, +0.11102230246251565E-15, +0.77715611723760958E-15, +0.44408920985006262E-15, +0.0000000000000000, +0.33306690738754696E-15, +0.11102230246251565E-15, +0.22204460492503131E-15, +0.55511151231257827E-15, +0.66613381477509392E-15, +0.66613381477509392E-15, +0.33306690738754696E-15, +0.33306690738754696E-15, +0.33306690738754696E-15, +0.33306690738754696E-15, +0.77715611723760958E-15, +0.33306690738754696E-15, +0.55511151231257827E-15, +0.55511151231257827E-15, +0.88817841970012523E-15
408reltol = sqrt(epsilon(1._TKC))
409reltol
410+0.14901161193847656E-7
411isClose(data, atad(1:size(data)), reltol = reltol)
412T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T
413call setAsserted(all(isClose(data, atad(1:size(data)), reltol = reltol)))
414
415lenData = getUnifRand(5, 100)
416lenData
417+44
418data = 1._TKC + getUnifRand(0._TKC, 1._TKC, lenData)
419data
420+1.34938812672340384579401274147850791, +1.04129816744660027917753004079894992, +1.80252115614277422018379647109665993, +1.42089486018824610048692675418010783, +1.30355310027292514045536653583353515, +1.38821554211189474153258611475470729, +1.06446632941455760614715737556925511, +1.11987785292360420895517182111211379, +1.47904611466330519327549166537057744, +1.19561222018827207433256875028358347, +1.59970065292563322209471115338551107, +1.12519589349938326146765524218378949, +1.36645195519043489606313831205602794, +1.87361330721245992848823248757978467, +1.79362435172835412562953623137883942, +1.03791371135703605581343155606469851, +1.56418683401054756568047356942139324, +1.52068801138771640404496998977946784, +1.81203701616421352903846856596022860, +1.41722665862368334133497144776866853, +1.40876993843963874059436715626257599, +1.10552740439643098884509910899615650, +1.15027577591546646903920473611474316, +1.70322255176954954689369435325148655, +1.74376837077953174415068275583713699, +1.06300677513864525514090492854493157, +1.80159093395531999345421881900864536, +1.96819713622665887907009243234575512, +1.24094499813886520799030414714581755, +1.98373587613086015017315800689415700, +1.40744618557630538017150262740801231, +1.64837235559526722583409514943219653, +1.32729068690252664269416275158834481, +1.72955199869148116549199170772838747, +1.30631266704282962166501961678227288, +1.22105819344195329691572493867941034, +1.73366549148559237275671794387354314, +1.39327842058785529539846566086110865, +1.54503549858153860513364597328615922, +1.88962705328568501838365920546104117, +1.84729495433125560132123253310840544, +1.66584933165749924537107642825209236, +1.66697021525675519004862274879742787, +1.24232234476154844684950532614521844
421
422atad = getFFTR(getFFTF(data))
423atad = atad * 2 / size(atad)
424atad
425+1.34938812672340384579401274147850791, +1.04129816744660027917753004079895012, +1.80252115614277422018379647109665955, +1.42089486018824610048692675418010783, +1.30355310027292514045536653583353553, +1.38821554211189474153258611475470749, +1.06446632941455760614715737556925511, +1.11987785292360420895517182111211379, +1.47904611466330519327549166537057744, +1.19561222018827207433256875028358347, +1.59970065292563322209471115338551088, +1.12519589349938326146765524218378911, +1.36645195519043489606313831205602813, +1.87361330721245992848823248757978487, +1.79362435172835412562953623137883942, +1.03791371135703605581343155606469851, +1.56418683401054756568047356942139324, +1.52068801138771640404496998977946804, +1.81203701616421352903846856596022802, +1.41722665862368334133497144776866853, +1.40876993843963874059436715626257599, +1.10552740439643098884509910899615650, +1.15027577591546646903920473611474296, +1.70322255176954954689369435325148655, +1.74376837077953174415068275583713661, +1.06300677513864525514090492854493176, +1.80159093395531999345421881900864516, +1.96819713622665887907009243234575512, +1.24094499813886520799030414714581755, +1.98373587613086015017315800689415739, +1.40744618557630538017150262740801231, +1.64837235559526722583409514943219653, +1.32729068690252664269416275158834520, +1.72955199869148116549199170772838747, +1.30631266704282962166501961678227288, +1.22105819344195329691572493867941054, +1.73366549148559237275671794387354352, +1.39327842058785529539846566086110865, +1.54503549858153860513364597328615922, +1.88962705328568501838365920546104117, +1.84729495433125560132123253310840564, +1.66584933165749924537107642825209236, +1.66697021525675519004862274879742806, +1.24232234476154844684950532614521844, -0.192592994438723585305597794258492732E-33, -0.288889491658085377958396691387739098E-33, +0.962964972193617926527988971292463659E-34, -0.288889491658085377958396691387739098E-33, +0.385185988877447170611195588516985464E-33, +0.962964972193617926527988971292463659E-34, -0.962964972193617926527988971292463659E-34, +0.962964972193617926527988971292463659E-34, -0.962964972193617926527988971292463659E-34, -0.192592994438723585305597794258492732E-33, +0.00000000000000000000000000000000000, -0.962964972193617926527988971292463659E-34, -0.192592994438723585305597794258492732E-33, +0.962964972193617926527988971292463659E-34, -0.288889491658085377958396691387739098E-33, +0.385185988877447170611195588516985464E-33, +0.288889491658085377958396691387739098E-33, -0.192592994438723585305597794258492732E-33, +0.288889491658085377958396691387739098E-33, +0.962964972193617926527988971292463659E-34
426reltol = sqrt(epsilon(1._TKC))
427reltol
428+0.138777878078144567552953958511352539E-16
429isClose(data, atad(1:size(data)), reltol = reltol)
430T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T
431call setAsserted(all(isClose(data, atad(1:size(data)), reltol = reltol)))
432
433lenData = getUnifRand(5, 100)
434lenData
435+10
436data = 1._TKC + getUnifRand((0._TKC, 0._TKC), (1._TKC, 1._TKC), lenData)
437data
438(+1.32873142, +0.439862311), (+1.20111871, +0.840382338), (+1.56901646, +0.281844258), (+1.64273715, +0.410495698), (+1.47816706, +0.730382621), (+1.63878024, +0.832411945), (+1.45539165, +0.693777204), (+1.65732288, +0.407264054), (+1.19527113, +0.388794065), (+1.33059156, +0.252408147)
439
440atad = getFFTR(getFFTF(data))
441atad = atad / size(atad)
442atad
443(+1.32873142, +0.439862341), (+1.20111871, +0.840382338), (+1.56901646, +0.281844288), (+1.64273715, +0.410495788), (+1.47816706, +0.730382442), (+1.63878024, +0.832411766), (+1.45539165, +0.693776846), (+1.65732288, +0.407263845), (+1.19527113, +0.388793975), (+1.33059144, +0.252408057), (+0.178813934E-6, -0.298023224E-7), (+0.119209290E-6, -0.149011612E-6), (-0.596046448E-7, +0.894069672E-7), (+0.00000000, +0.596046448E-7), (+0.596046448E-7, +0.268220901E-6), (+0.596046448E-7, +0.149011612E-6)
444reltol = sqrt(epsilon(1._TKC))
445reltol
446+0.345266977E-3
447isClose(data, atad(1:size(data)), reltol = reltol)
448T, T, T, T, T, T, T, T, T, T
449call setAsserted(all(isClose(data, atad(1:size(data)), reltol = reltol)))
450
451lenData = getUnifRand(5, 100)
452lenData
453+25
454data = 1._TKC + getUnifRand((0._TKC, 0._TKC), (1._TKC, 1._TKC), lenData)
455data
456(+1.6853137614995739, +0.41385034054789194), (+1.7192283296887838, +0.90666468046146664), (+1.7127750558941583, +0.57440698159126780), (+1.8219746646370685, +0.24649816212012132), (+1.3780797656873469, +0.79919291241394075), (+1.3007620085189222, +0.96880729510480101), (+1.5786123958166942, +0.35747898901812092), (+1.5055786673718430, +0.75780299603212586), (+1.0453624456914787, +0.65318406508616300), (+1.7417427316017200, +0.26672604810570044), (+1.7861973387644967, +0.49633051671795525), (+1.3339354595846187, +0.86594552477473807), (+1.0757725216696086, +0.69851622584928696), (+1.4843387831293833, +0.28345664402316095), (+1.8315738499540468, +0.46586990650651083E-2), (+1.0091880851603010, +0.65251446718728046), (+1.8429401292214354, +0.55398134206304916), (+1.9901248163208285, +0.63186064164864608), (+1.3550104034392021, +0.35146544415009617), (+1.2911042483750086, +0.67179230166192772), (+1.2959823691637582, +0.21666721503528907), (+1.5255825696305432, +0.20013618755144391), (+1.6621977024708974, +0.14937889341865884), (+1.7954192676868641, +0.22775897427078107), (+1.3388400232080480, +0.28436791961006525)
457
458atad = getFFTR(getFFTF(data))
459atad = atad / size(atad)
460atad
461(+1.6853137614995739, +0.41385034054789205), (+1.7192283296887838, +0.90666468046146687), (+1.7127750558941586, +0.57440698159126768), (+1.8219746646370685, +0.24649816212012157), (+1.3780797656873469, +0.79919291241394075), (+1.3007620085189218, +0.96880729510480090), (+1.5786123958166942, +0.35747898901812108), (+1.5055786673718425, +0.75780299603212564), (+1.0453624456914787, +0.65318406508616311), (+1.7417427316017196, +0.26672604810570061), (+1.7861973387644967, +0.49633051671795497), (+1.3339354595846185, +0.86594552477473785), (+1.0757725216696086, +0.69851622584928674), (+1.4843387831293828, +0.28345664402316095), (+1.8315738499540464, +0.46586990650650806E-2), (+1.0091880851603010, +0.65251446718728023), (+1.8429401292214354, +0.55398134206304916), (+1.9901248163208287, +0.63186064164864586), (+1.3550104034392023, +0.35146544415009628), (+1.2911042483750088, +0.67179230166192749), (+1.2959823691637582, +0.21666721503528907), (+1.5255825696305432, +0.20013618755144369), (+1.6621977024708974, +0.14937889341865868), (+1.7954192676868641, +0.22775897427078096), (+1.3388400232080477, +0.28436791961006530), (+0.33306690738754696E-15, +0.55511151231257827E-16), (+0.66613381477509392E-15, +0.55511151231257827E-16), (+0.44408920985006262E-15, +0.0000000000000000), (+0.44408920985006262E-15, +0.27755575615628914E-15), (+0.55511151231257827E-15, +0.44408920985006262E-15), (+0.55511151231257827E-15, +0.19428902930940239E-15), (+0.55511151231257827E-15, +0.61062266354383610E-15)
462reltol = sqrt(epsilon(1._TKC))
463reltol
464+0.14901161193847656E-7
465isClose(data, atad(1:size(data)), reltol = reltol)
466T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T
467call setAsserted(all(isClose(data, atad(1:size(data)), reltol = reltol)))
468
469lenData = getUnifRand(5, 100)
470lenData
471+23
472data = 1._TKC + getUnifRand((0._TKC, 0._TKC), (1._TKC, 1._TKC), lenData)
473data
474(+1.70253227264203933354580566144284279, +0.916068137873386138340332510277841468), (+1.09978148417094498689678724178837807, +0.185976538384103642109087159579709017), (+1.47680145450847586516040786424969544, +0.745602462288547680046281309253149569), (+1.69734547741612571869218280755041424, +0.801132103773256744551060498502075084), (+1.27065105889399158923816404279451732, +0.810195930865451614129686684402284348), (+1.07139550837028317914572604009754028, +0.378733653747373202155746492181216939), (+1.74665217477180743011625134758389928, +0.197268301291169599828232766723199875E-1), (+1.37023642850788612445280257304447263, +0.736623489287089318622461458438794760), (+1.06092015009561205953465626666221046, +0.249942238028434693024783175483303420), (+1.47438319685962579635460474209334740, +0.717854833387468334471259445580442603), (+1.54007931379121673442623765874316359, +0.154120004643801810539674316824232428), (+1.85318153761539956745515630130107018, +0.378914395057125941986787134837651732), (+1.73598358694466476783793755817779827, +0.474304500528082994364132212991559888), (+1.95724866282515816925237066378702481, +0.805641333675432459161439404900351050), (+1.19040068937059092673387506154802630, +0.742254830761816605626688053546280726), (+1.10866230215622386093520433255781061, +0.824269872808710046594987487847358893), (+1.86652428917290774881493516749223584, +0.506371190438013460121682391425904716), (+1.41914666495768334462559386758294986, +0.295089163758531501451119211801190914), (+1.89453936998378403670206484358721501, +0.785439209590517637958718821540831848), (+1.02534126177448910034653504797705129, +0.742435376400872948095452628032104432), (+1.83844721649532100814443504041215611, +0.624622267356563421556192541130271251), (+1.38007313793134957851292400572032912, +0.472591299152884234427934348801635214), (+1.20188944878971843633886696311817031, +0.446781267044035858891059475054489937)
475
476atad = getFFTR(getFFTF(data))
477atad = atad / size(atad)
478atad
479(+1.70253227264203933354580566144284260, +0.916068137873386138340332510277841468), (+1.09978148417094498689678724178837827, +0.185976538384103642109087159579709065), (+1.47680145450847586516040786424969564, +0.745602462288547680046281309253149473), (+1.69734547741612571869218280755041424, +0.801132103773256744551060498502074988), (+1.27065105889399158923816404279451732, +0.810195930865451614129686684402284348), (+1.07139550837028317914572604009754067, +0.378733653747373202155746492181216842), (+1.74665217477180743011625134758390005, +0.197268301291169599828232766723197227E-1), (+1.37023642850788612445280257304447282, +0.736623489287089318622461458438794471), (+1.06092015009561205953465626666221027, +0.249942238028434693024783175483303420), (+1.47438319685962579635460474209334759, +0.717854833387468334471259445580442410), (+1.54007931379121673442623765874316359, +0.154120004643801810539674316824232332), (+1.85318153761539956745515630130107018, +0.378914395057125941986787134837651732), (+1.73598358694466476783793755817779807, +0.474304500528082994364132212991559792), (+1.95724866282515816925237066378702481, +0.805641333675432459161439404900350954), (+1.19040068937059092673387506154802630, +0.742254830761816605626688053546280823), (+1.10866230215622386093520433255781061, +0.824269872808710046594987487847359086), (+1.86652428917290774881493516749223564, +0.506371190438013460121682391425904620), (+1.41914666495768334462559386758294986, +0.295089163758531501451119211801190818), (+1.89453936998378403670206484358721540, +0.785439209590517637958718821540832040), (+1.02534126177448910034653504797705129, +0.742435376400872948095452628032104625), (+1.83844721649532100814443504041215630, +0.624622267356563421556192541130271251), (+1.38007313793134957851292400572032950, +0.472591299152884234427934348801635117), (+1.20188944878971843633886696311817050, +0.446781267044035858891059475054489840), (+0.00000000000000000000000000000000000, +0.144444745829042688979198345693869549E-33), (+0.00000000000000000000000000000000000, -0.962964972193617926527988971292463659E-34), (+0.577778983316170755916793382775478196E-33, +0.481482486096808963263994485646231830E-34), (-0.962964972193617926527988971292463659E-34, -0.192592994438723585305597794258492732E-33), (+0.962964972193617926527988971292463659E-34, -0.288889491658085377958396691387739098E-33), (+0.00000000000000000000000000000000000, +0.00000000000000000000000000000000000), (+0.385185988877447170611195588516985464E-33, +0.192592994438723585305597794258492732E-33), (-0.962964972193617926527988971292463659E-34, -0.192592994438723585305597794258492732E-33), (+0.288889491658085377958396691387739098E-33, -0.288889491658085377958396691387739098E-33)
480reltol = sqrt(epsilon(1._TKC))
481reltol
482+0.138777878078144567552953958511352539E-16
483isClose(data, atad(1:size(data)), reltol = reltol)
484T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T
485call setAsserted(all(isClose(data, atad(1:size(data)), reltol = reltol)))
486
487lenData = getUnifRand(5, 100)
488lenData
489+48
490data = 1._TKC + getUnifRand(0._TKC, 1._TKC, lenData)
491data
492+1.47820723, +1.11837649, +1.11516261, +1.95121765, +1.91581559, +1.90463805, +1.04372644, +1.55271912, +1.25922394, +1.36798668, +1.22661173, +1.80899692, +1.05353165, +1.89220643, +1.33239365, +1.10004997, +1.64970863, +1.07824540, +1.91912127, +1.34425116, +1.71818483, +1.04063261, +1.71735489, +1.12263942, +1.36682153, +1.17146397, +1.47653604, +1.62502265, +1.25026250, +1.18942523, +1.30306244, +1.47393548, +1.68877149, +1.46353841, +1.76702213, +1.44114530, +1.07454109, +1.67786455, +1.17278242, +1.26247239, +1.47087479, +1.40837371, +1.46905053, +1.75410175, +1.73818159, +1.39765787, +1.49998093, +1.15985107
493
494atad = getFFTR(getFFTF(data))
495atad = atad * 2 / size(atad)
496atad
497+1.47820711, +1.11837637, +1.11516261, +1.95121765, +1.91581535, +1.90463805, +1.04372644, +1.55271912, +1.25922382, +1.36798668, +1.22661173, +1.80899680, +1.05353165, +1.89220631, +1.33239365, +1.10004985, +1.64970863, +1.07824540, +1.91912115, +1.34425116, +1.71818483, +1.04063272, +1.71735477, +1.12263942, +1.36682129, +1.17146397, +1.47653580, +1.62502253, +1.25026250, +1.18942523, +1.30306220, +1.47393537, +1.68877149, +1.46353829, +1.76702189, +1.44114542, +1.07454085, +1.67786431, +1.17278242, +1.26247215, +1.47087467, +1.40837359, +1.46905077, +1.75410163, +1.73818159, +1.39765775, +1.49998093, +1.15985119, -0.119209290E-6, +0.00000000, +0.119209290E-6, +0.417232513E-6, -0.119209290E-6, +0.00000000, +0.178813934E-6, +0.178813934E-6, -0.178813934E-6, +0.238418579E-6, +0.596046448E-7, +0.357627869E-6, +0.238418579E-6, +0.238418579E-6, +0.178813934E-6, +0.298023224E-6
498reltol = sqrt(epsilon(1._TKC))
499reltol
500+0.345266977E-3
501isClose(data, atad(1:size(data)), reltol = reltol)
502T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T
503call setAsserted(all(isClose(data, atad(1:size(data)), reltol = reltol)))
504
505lenData = getUnifRand(5, 100)
506lenData
507+52
508data = 1._TKC + getUnifRand(0._TKC, 1._TKC, lenData)
509data
510+1.4292635796862330, +1.6867779814888100, +1.4590659056816231, +1.5512800096950263, +1.6179283418710786, +1.0929363619488970, +1.5770965277333344, +1.7981749889109795, +1.3160442690232776, +1.5632116592537078, +1.0563879890286039, +1.9437201814771941, +1.9345605262332495, +1.5628603613350442, +1.2436095162940566, +1.0126356891765815, +1.8368294794243989, +1.9062187012832088, +1.6529000195754469, +1.6734781392633733, +1.0494417878477256, +1.9084377940764221, +1.7971676400780223, +1.7694109188218132, +1.1707480188561568, +1.5877869975744132, +1.0559950445470876, +1.4680121653305878, +1.3537046380886073, +1.4244124863095202, +1.9365895039434626, +1.4425672477721991, +1.3940105910846696, +1.8162178969987341, +1.2991614579043256, +1.0409114814717462, +1.6035737194017554, +1.2927363789846020, +1.1223963881901988, +1.6018284898459214, +1.5157140522560346, +1.6195890091009379, +1.4311902944347579, +1.9063052733677792, +1.6342609904572347, +1.4934755872638066, +1.5131030897595332, +1.2640729142372602, +1.9797067827602832, +1.9231622952701888, +1.0743022975296930, +1.1162395833125607
511
512atad = getFFTR(getFFTF(data))
513atad = atad * 2 / size(atad)
514atad
515+1.4292635796862327, +1.6867779814888102, +1.4590659056816229, +1.5512800096950270, +1.6179283418710786, +1.0929363619488970, +1.5770965277333340, +1.7981749889109793, +1.3160442690232772, +1.5632116592537073, +1.0563879890286039, +1.9437201814771938, +1.9345605262332488, +1.5628603613350442, +1.2436095162940564, +1.0126356891765811, +1.8368294794243987, +1.9062187012832086, +1.6529000195754462, +1.6734781392633733, +1.0494417878477253, +1.9084377940764217, +1.7971676400780219, +1.7694109188218126, +1.1707480188561563, +1.5877869975744126, +1.0559950445470874, +1.4680121653305873, +1.3537046380886069, +1.4244124863095200, +1.9365895039434615, +1.4425672477721982, +1.3940105910846696, +1.8162178969987341, +1.2991614579043254, +1.0409114814717462, +1.6035737194017554, +1.2927363789846020, +1.1223963881901988, +1.6018284898459210, +1.5157140522560342, +1.6195890091009375, +1.4311902944347579, +1.9063052733677790, +1.6342609904572343, +1.4934755872638066, +1.5131030897595330, +1.2640729142372602, +1.9797067827602830, +1.9231622952701890, +1.0743022975296930, +1.1162395833125609, +0.22204460492503131E-15, -0.11102230246251565E-15, +0.66613381477509392E-15, +0.0000000000000000, +0.0000000000000000, +0.0000000000000000, +0.22204460492503131E-15, +0.55511151231257827E-15, +0.11102230246251565E-15, +0.33306690738754696E-15, +0.66613381477509392E-15, +0.88817841970012523E-15
516reltol = sqrt(epsilon(1._TKC))
517reltol
518+0.14901161193847656E-7
519isClose(data, atad(1:size(data)), reltol = reltol)
520T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T, T
521call setAsserted(all(isClose(data, atad(1:size(data)), reltol = reltol)))
522
523lenData = getUnifRand(5, 100)
524lenData
525+7
526data = 1._TKC + getUnifRand(0._TKC, 1._TKC, lenData)
527data
528+1.86833938479387488060918321198483646, +1.24560079322233101466777215106078798, +1.20389353800237589796619295555184076, +1.02292359311748096752353194371658972, +1.17447163718773691329888164238611098, +1.82717189178758413457985022063151538, +1.75212615430796006938015817951666628
529
530atad = getFFTR(getFFTF(data))
531atad = atad * 2 / size(atad)
532atad
533+1.86833938479387488060918321198483646, +1.24560079322233101466777215106078778, +1.20389353800237589796619295555184076, +1.02292359311748096752353194371658972, +1.17447163718773691329888164238611117, +1.82717189178758413457985022063151538, +1.75212615430796006938015817951666628, -0.962964972193617926527988971292463659E-34
534reltol = sqrt(epsilon(1._TKC))
535reltol
536+0.138777878078144567552953958511352539E-16
537isClose(data, atad(1:size(data)), reltol = reltol)
538T, T, T, T, T, T, T
539call setAsserted(all(isClose(data, atad(1:size(data)), reltol = reltol)))
540
541
Test:
test_pm_fftnr


Final Remarks


If you believe this algorithm or its documentation can be improved, we appreciate your contribution and help to edit this page's documentation and source file on GitHub.
For details on the naming abbreviations, see this page.
For details on the naming conventions, see this page.
This software is distributed under the MIT license with additional terms outlined below.

  1. If you use any parts or concepts from this library to any extent, you must acknowledge the usage by citing the relevant publications of the ParaMonte library.
  2. If you regenerate any parts/ideas from this library in a programming environment other than those currently supported by this ParaMonte library (i.e., other than C, C++, Fortran, MATLAB, Python, R), you must also ask the end users to cite this original ParaMonte library.
Author:
Fatemeh Bagheri, Tuesday 11:34 PM, August 10, 2021, Dallas, TX

Definition at line 611 of file pm_fftnr.F90.


The documentation for this interface was generated from the following file: